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13  ABSTRACT 


The  objective  of  this  program  was  to  develop  a  computer  program  capable  of 
accurately  predicting  the  internal  ballistics  of  motors  designed  for  operation 
without  nozzles.  To  accomplish  this,  previously  obtained  data  were  examined 
in  detail.  In  addition,  12  motors  were  fabricated  and  test  fired  to  augment  the 
range  of  available  test  data.  The  computer  program  developed  had  to  consider 
the  effects  of  erosive  burning,  grain  deflection,  and  incomplete  metal  combustion 
upon  nozzleless  motor  performance.  Input  requires  knowledge  of  motor  geometry, 
propellant  thermochemistry,  base  burn  rate,  and  specification  of  erosive  burning 
and  grain  deflection  parameters.  Comparisons  of  predicted  and  actual  results 
show  that  where  an  adequate  data  base  exists,  with  regard  to  erosive  burning  and 
grain  deflection  effects,  accurate  predictions  will  be  provided  by  the  computer 
program.  This  applies  to  all  performance  parameters,  including  instantaneous 
thrust  and  pressure,  burn  time,  and  specific  and  total  impulse.  Even  lacking 
such  a  data  base,  the  program  will  correctly  predict  the  overall  parameters  of 
specific  impulse  and  total  impulse.  To  fully  utilize  the  capabilities  provided  by 
the  computer  program,  further  study  of  erosive  burning  and  grain  deflection  is 
required.  The  technical  program  is  described  in  AFRPL-TR-7  3- i9. 
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FOREWORD 


This  technical  report  summarizes  work  performed  from  1  June  1972  to 
31  January  1973  under  Contract  No.  F0461 1  -72-C-0076  by  Lockheed 
Propulsion  Company  (LPC),  Redlands,  California,  for  the  Air  Force  Rocket 
Propulsion  Laboratory  (AFRPL),  Edwards,  California. 

The  work  reported  herein  was  performed  unoer  the  technical  direction 
of  Lt.  H.  Barbarika  of  the  Air  Force  Rocket  Propulsion  Laboratory.  The 
Lockheed  Propulsion  Company  program  manager  was  Mr.  D.  E.  Taylor, 
and  the  project  engineer  was  Mr.  K.  R.  Small.  Mr.  D.  P.  Harry  conducted 
computer  program  development  efforts,  and  Mr.  C.  F.  Price  interpreted 
ballistic  results  and  cinefluorographic  test  data.  Propellant  characterization 
and  processing  were  under  the  direction  of  Mr.  I.  L.  Markovitch  and  Mr. 

F.  C.  Anderson.  The  Lockheed  Propulsion  Company  project  number  was 
MPO  676.  / 

This  technical  report  has  been  reviewed  and  is  approved. 


Charles  R.  Cooke 

Chief,  Solid  Rocket  Division 

AFRPL/MK 
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SECTION  I 

INTRODUCTION  AND  SUMMARY 


The  potential  of  a  solid  rocket  operating  without  an  exit  nozzle  for 
low-cost  or  high-altitude  applications  has  been  known  for  some  time. 

Recently,  the  feasibility  of  the  concept  was  established  in  a  test  program 
conducted  by  United  Technology  Center  (UTC)  on  Contract  F046 1 1 -7  1 -C -0050. 
Reproducibility  of  performance  was  demonstrated  and  grain  shape  variations 
that  delivered  increased  performance  were  identified.  The  data  acquired 
provided  a  good  base  for  establishment  of  an  analytical  model  to  accurately 
predict  the  ballistic  performance  of  nozzleless  motors.  Such  a  model  would 
provide  methodology  to  maximize  the  potential  of  such  a  rocket  by  design, 
and  to  permit  evaluation  of  modifications  that  will  optimize  performance  at 
no  sacrifice  in  simplicity  of  fabrication. 

The  foundation  for  the  computer  program  was  a  nozzleless  rocket 
ballistics  computer  program  already  existing  at  LPC,  which  had  evolved  in 
sophistication  over  the  year  preceding  this  contract.  The  program  had  its 
inception  in  a  conventional  ballistics  program  with  four  descriptions  of  ero¬ 
sive  burning  included,  and  with  allowances  to  treat  axial  variations  with 
steady-state  versions  of  the  one -dimensional  channel  flow  equations.  Included 
in  the  original  ballistics  program  was  a  projection  of  the  grain  profile  shape. 
The  program  was  modified  to  make  it  applicable  to  the  fast  transients  that 
occur  in  ignition  and  tailoff.  For  chr  purpose  of  examining  ignition  and  tailoff 
transients,  the  usual  lumped  volume  approaches,  which  have  been  developed 
to  apply  to  conventional  rockets,  are  not  applicable  to  the  nozzleless  rocket. 
This  is  because  they  are  incapable  of  treating  the  impr  rtant  axial  variations 
in  the  nozzleless  rocket,  which  dominate  its  behavior.  The  program  modifi¬ 
cation  was  developed  from  an  examination  of  the  pertinent,  unsteady,  one¬ 
dimensional  flow  equations.  It  was  found  that  an  approximate  method  of 
solution  to  the  transient  problem  could  be  readily  obtained.  Briefly,  it 
retains  a  time-dependent  term  in  the  continuity  equacion  only  while  retaining 
all  important  axially  dependent  terms. 

Computer  program  input  requires  detailed  knowledge  of  motor 
geometry  features  and  propellant  thermochemistry.  In  addition,  the  user 
must  select  inputs  related  to  erosive  burning,  grain  deformation,  and  metal 
combustion.  The  output  provides  all  ove  all  performance  parameters  (head- 
end  pressure,  thrust,  mass  flow  rate,  etc),  as  well  as  a  description  of 
internal  motor  geometry,  all  versus  time.  Additional  output  is  available  as 
a  user  option.  Serial  runs  only  require  input  of  changes.  Computer  run 
time  varies  between  10  seconds  and  1  minute  on  a  CDC  6400  machine. 

This  User's  Manual  completely  describes  and  documents  the 
nozzleless  rocket  motor  internal  ballistics  computer  program.  Test  cases 
arc  included  to  demonstrate  the  use  of  the  program.  The  model  is 
operational  on  the  AFRPL  computer  facility. 
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SECTION  II 
ANALYSIS 


1.  FLOW  EQUATIONS 

The  equations  for  one-dimensional  flow  used  in  the  computer  pro¬ 
gram  are  based  on  the  following  assumptions: 

•  Particles  occupy  negligible  volume  compared  to  the  gas 

•  Particle  velocity  is  a  constant  fraction  of  gas  velocity 

•  The  gas  follows  perfect  gas  relationships 

•  Mass  added  from  the  propellant  surface  has  no  axial 
component  of  velocity 


With  these  assumptions,  the  equations  that  apply  are  (symbols  are  defined 
in  Appendix  A) : 

STATE 

p  =  p(l-  <  )Rg  T 

CONTINUITY 

£<PA>  +  &  H'-O-O  1  UBA]  = 

MOMENTUM 

PA  [l-.(l-f«)J  Ug  + 

Ae  Ug  [  *-(>-0  J 

+  A 

-|£-  +  F  =0 

OX  X 

ENERGY 
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i 

Choking  occurs  when  the  gas  phase 
the  apparent  sonic  velocity,  c,  given  by 

velocity,  U  .  becomes  equal  to 

c  =  Ju  l-<)  Rg  T 

C 
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where: 

c  =  specific  heat  at  constant  pressure  of  the  mixture 
p  (for  f  f  1) 

c  =  specific  heat  at  constant  volume  of  the  mixture 
(for  f  j.  1) 

The  mass  fraction  of  particles  in  the  flow  is  higher  wh.?i  the 

particles  lag  behind  the  gas,  and  can  be  described  in  terms  of  th:  amount 

of  lag  and  t  : 

°  o 

t 

<  =  -  ° - 

f  +  d-f)  •„ 

2.  EROSIVE  BURNING  CORRELATIONS 

The  computer  program  presently  incorporates  an  erosive  burn 
rate  description  of  the  form 

r  -  r„  +  a  <G”  -  00")  (af 

where: 

m  =  0.8 

c  =  sonic  velocity 

a,  Gq  =  constants  determined  by  the  propellant 

To  provide  the  user  with  some  generality,  the  program  also 
retains  the  Lenoir-Robillard  form  of  erosive  burn  rate  description 


r 

where: 


_  .  a  (G)0' 

ro  +  -ib&r 


'-¥) 


p  =  solid  propellant  density 
s 

a,  P  =  constants  determined  by  the  propellant  (a  related 
theoretically  to  a  heat  balance) 
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3.  ALUMINUM  COMBUSTION  AND  SLIP  FLOW 

Empirical  equations  are  used  to  fit  results  of  individual  particle 
calculations.  The  correlations  used  are: 


F  = 


\0. 2271 


u 


A1 


A1  “  V, 


U, 


=  ZAhQi 

A12Q3  ~  Vf 


/  t  \0,6633  /  v  0 . 1 1  51 

{■k)  (sfe)  U-0.0084  Do) 

0.5534fe)-(^)° 

(mJ 


V 0.0842 


\  0.02 34 


=  0.9470 


where: 


/  t)  \  0 .  0  38 1 

UooJ 


L  =  motor  length,  in. 

p  =  head-end  pressure,  psia 

D  =  initial  diameter  of  aluminum  agglomerate  leaving 
propellant  surface,  microns 

F  =  fraction  of  aluminum  burned 

U^i  =  velocity  of  aluminum  particle  at  motor  exit,  ft/sec 


U 


fA1  =  velocity  of  aluminum  oxide  particle  at  motor  exit,  ft/sec 

Ai2U3 


Vf  =  gas  velocity  at  motor  exit,  ft/sec 


4.  THERMOCHEMISTRY  INTERPRETATION 


Thermochemical  properties  of  the  combustion  products  determined 
from  standard  equilibrium  thermochemistry  programs  are  assumed  input 
to  the  model  (the  input  form  is  described  in  Section  IV).  The  input  should 
cover  at  least  the  range  of  the  fraction  of  aluminum  burned  expected  in  the 
particular  motor,  and  is  obtained  from  so-called  T*  runs  of  the  thermo¬ 
chemistry  program,  meaning  that  a  fraction  of  the  aluminum  is  treated  as 
a  "non-reacting  specie".  The  theoretical  thermochemistry  program  in  use 
at  AFRPL  may  be  obtained  from 

Robert  A.  Hall  (714)  553-2196 

AFRPL/ TETT 

Edwards  Air  Force  Base,  CA  93523 
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Calculation  of  gas  properties  from  thermochemical  program  inputs, 
the  fraction  of  aluminum  burned,  and  slip  flow  correlations  employs  the 
following  relationships: 

Rg  =  Ru/Mg 

«  =  1  -  Mg  (Moles  gas)/lOO 

...  +  fAl,Q,> 

-  *ai2o,  +  <o 


where  A1  is  the  fraction  of  aluminum  in  the  propellant. 
V  =  Cp,  p/(Cp,  p  -  Rg  M°}“  ) 


€  = 


f  l1  -«o)  +  *o 
Rp  =  (l  - «)  Rg 


v«  "  Cp' 


C £ 


Cp 

5.  GRAIN  DEFLECTION 


-  1.987 

_  Moles  prod 
=  Cp.g  +  (Cp,  p  MoIes  £as 


-  Cp,g)/f 


Grain  deflection  is  considered  in  the  computer  program  as  a  product 
of  steady-state  deflection  modified  by  the  ratio  of  transient  to  steady-state 
motion.  The  steady -state  deflection  is  obtained  as  a  function  of  motor  length 
from  computer  solutions  to  the  boundary  value  problem  for  the  specific  pro¬ 
pellant  geometry  of  interest.  The  results  are  linear  with  head  pressure  to 
ambient  pressure  difference  for  a  given  pressure  distribution  in  the  port. 
Grain  deflection  has  been  found  to  be  nearly  linear  with  the  web  fraction 
burned  for  a  given  AP,  ana  this  approximation  is  used  in  the  model: 


Ar 


head 


ref 


v/eb 
web,  l 


where  the  steady- state  deflection  for  the  initial  web  and  a  given  reference 
pressure  is  input  as  a  function  of  length.  The  throat  web  fraction  is  used  as 
representative  of  the  important  effects  on  motor  operation. 
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The  ra  .o  of  transient  to  steady-state  deflection  is  determined  by 
assuming  the  grain  to  act  as  a  critically  damped  spring-mass  system  driven 
by  a  pressure  ramp: 

x  +  2c>  x  +  u  2  x  =  Pt/m  +  P./m 
n  n  '0' 

where 


n  m 

The  integrated  equation  for  an  increment  of  calculation  where 
initial  conditions  are  not  zero  is 


K  (1  -  unt)  xQ  +  (2un  x0  +  x0)  Kt  -  P0  (1  +  wnt) 

+  P(2/u,n+t)  +  P0  +  P  (t  -  2/%) 

and  the  initial  rate  x  is  determined  by  differentiating  x(t).  The  deflection 
equation  is  integrated  with  respect  to  time  in  the  model  in  closed  form  for 
an  increment  t  =  At  but  is  iterated  as  P  is  iterated.  Values  of  the  input 
coefficients  K  and  are  selected  to  fit  experimental  data. 

6.  CHOKING  CONSTRAINTS 

Iterative  calculation  to  satisfy  choking  constraints  is  by  far  the 
most  time-consuming  operation  within  the  nozzleless  motor  ballistic  pre¬ 
diction  program.  Consequently,  sophisticated  correction  terms  are  gene¬ 
rated  to  predict  successive  trial  values.  The  logic  used  for  the  special 
cases  that  apply  to  nozzleless  motors  is  summarized  in  the  following 
paragraphs. 


*(t)  _ 
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head 
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a.  Choked  Flow,  Throat  at  Port  Exit 

A  head-end  pressure  P(l)  is  selected  and  the  integration  down 
the  port  is  performed.  If  the  velocity  of  the  throat  is  less  than  the  local 
sonic  velocity,  an  adjustment  is  made  to  decrease  head-end  pressure,  and 
the  iteration  continues.  If,  on  the  other  hand,  the  velocity  exceeds  the  sonic 
velocity,  a  trap  is  executed  that  limits  velocity  and  generates  an  error  signal 
by  "breaking"  the  momentum  equation.  An  adjustment  term  is  computed  to 
increase  P(l),  and  iteration  continues  until  the  specified  convergence  is 
obtained. 


b.  Subsonic  Flow 

If  the  entire  motor  is  operating  with  subsonic  flow,  the  head- 
end  pressure  is  adjusted  to  cause  exit  static  pressure  to  match  ambient 
pressure,  and  again,  an  approximate  relationship  is  used  to  calculate 
successive  trial  values. 
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c.  Choked  Flow,  Throat  m  Port 

Where  the  motor  contacts  an  exit  cone,  the  throat  is  located 
initially  near  the  nominal  minimum  area.  Solution  for  the  subsonic  flow 
portion  of  the  motor  up  to  the  throat  is  the  same  as  previously  described. 
The  iteration  is  completed  to  match  the  throat  constraint  and  then  the  super 
sonic  portion  is  calculated  using  the  same  equations  but  with  modified 
convergence  logic  because  the  computational  stability  is  altered.  Two 
added  sets  of  logic  are  checked: 

(i)  If  the  flow  expands  supersonically  to  the  exit,  the 
problem  is  complete.  The  possibility  that  exit  static 
pressure  is  sufficiently  lower  than  ambient  pressure 

to  allow  a  fully  expanded  nozzle  without  separation 
is  ignored  in  the  model.  The  throat  regression 
is  faster  than  that  in  the  exit  cone  and  the  expansion 
ratio  dec  reases ,  causing  (ii),  below,  but  not  an 
overexpanded  nozzle. 

(ii)  If  the  flow  is  calculated  io  decelerate  below  the  sonic 
velocity,  the  velocity  is  again  constrained  and  an 
error  signal  is  generated.  Note  that  the  mass  addi¬ 
tion  in  the  port  causes  a  second  "throat"  of  signifi¬ 
cantly  larger  area  than  the  first,  and  the  minimum 
area  may  operate  in  subsonic  flow.  A  sequential 
recovery  procedure  is  required.  The  "throat"  is 
first  identified  with  the  exit  of  the  port.  The  integra¬ 
tion  of  subsonic  flow  down  the  port  is  repeated  and 
the  new  throat  constraint  location  is  identified  by  the 
criteria: 

d  A  w 

_ £  .  at! tied  V  *  1  2  t 

A  '  w  1  ♦  <0  *  p 

P 

where  the  barred  quantities  are  the  average  for  the 
increment.  Thus,  the  "throat"  is  located  in  the 
expansion  region  at  the  node  where  the  area  increase 
just  overrides  the  effects  of  mass  addition  and  friction. 

A  note  in  program  output  is  printed  at  each  attempt  to 
relocate  the  throat.  The  tentative  solution  is  then 
iterated  to  the  specified  convergence  accuracy. 


11-6 


7 


AFRPL-TR- 73-20 


SECTION  III 
PROGRAM  STRUCTURE 


The  nozzleless  rocket  motor  ballistics  model  is  designed  to  compute 
a  time-history  of  motor  behavior  during  the  simulated  run.  The  logic  of 
solving  the  equations  presented  in  Section  II  is  described  in  this  section 
followed  by  a  discussion  of  the  computer  program  subroutine  structure.  The 
details  of  program  input  and  output  are  given  in  Sections  IV  and  V. 

1.  SOLUTION  LOGIC 

The  general  scheme  employed  in  solving  the  equations  of  one¬ 
dimensional  flow  is  shown  in  simplified  terms  in  Figure  III- I.  The  pro¬ 
cedure  that  follows  the  initialization  of  a  new  case  is  as  follows: 

•  A  head-end  pressure  is  assumed.  The  calculations  for  a 
single  axial  increment  are  made  and  iterated  to  completion. 

The  successive  increments  down  the  port  are  then  calculated 
until  the  choke  point  or  exit  is  reached. 

•  At  the  choke  point,  the  local  Mach  number  is  checked.  If 
the  calculated  value  is  not  within  tolerance,  a  correction 

to  head-end  pressure  is  made  and  the  calculation  is  restarted 
at  the  head.  As  successive  iterations  are  made  to  satisy  the 
choke  constraint  (outer  loop),  the  axial  stations  (inner  loop) 
are  continuously  iterated.  The  same  logic  applies  if  the  throat 
is  unchoked  except  that  p  .  =  p  ,  is  the  constraint. 

•  If  a  supersonic  exit  cone  exists  downstream  of  the  throat,  the 
axial  increments  are  continued  to  the  end  of  the  motor. 

•  Parameter'  used  in  transient  terms  and  integrating  calcula¬ 
tions  are  saved,  and  the  time  increment  is  advanced.  The 
calculations  resume  with  a  new  head-end  pressvire  assumption. 

•  As  the  calculations  reach  the  end  of  a  specified  run  time,  or 
the  motor  has  burned  out,  the  input  for  the  next  case  is  read 
and  the  new  case  is  started. 

""tc  computational  logic,  thus,  is  arranged  as  a  triple-loop  structure  with  the 
:ial  increment  (inner  loop),  choke  constraint  (outer  loop),  and  time  incre¬ 
ments  (stepping). 

A  more  detailed  diagram  of  program  logic  is  shown  in  Figure  III- ^ . 
This  representation  is  still  simplified  with  regard  to  the  individual  branching 
for  special  cases,  but  is  intended  to  describe  the  function  of  groupings  of 
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program  stops  internal  to  the  program,  and  presumes  valid  input.  Figure 
III-2  is  organized  in  the  same  order  as  the  source  program. 

•  The  compiled-in  constants  are  established  and  some  core  areas 
are  cleared  to  zero.  The  program  does  not  require  a  clear 
core . 

•  A  heading  card  is  read  and  printed.  An  integer  is  required 
in  CC1  to  indicate  an  ^er  case;  otherwise,  a  normal  exit  is 
interpreted. 

•  A  subroutine  to  control  the  reading  and  organizing  of  tabular 
data  is  called  by  the  main  program.  The  remaining  data, 
integers  and  real  data  input,  are  returned  in  intermediate 
storage  which  is  inspected  for  non-zero  (nz)  values.  Any 

nz  values  are  then  placed  in  working  storage  and  overwrite 
the  previous  values. 

•  Logical  interpretations  are  made  of  the  input  data;  for  example, 
the  throat  is  assumed  at  the  exit  unless  an  input  node  number 

is  given. 

•  Plenum  and  throat  geometry  for  a  nozzled  motor  is  calculated. 
These  steps  are  degenerate  unless  input  data  is  provided. 

•  Interpretation  of  one  or  two  propellant  cases  is  made  and  some 
constants  are  computed.  A  printed  note  identifies  the  pro¬ 
pellant  and  erosive  burn  model  to  be  used. 

•  The  geometry  of  the  grain  is  established  from  tabular  input  of 
inner  diameter,  outer  diameter,  and  axial  step  size.  Also, 
some  parameters  are  cleared  or  initialized  to  input  values. 

•  A  summary  of  the  key  words/data  that  describe  program 
options  is  printed.  For  example,  the  input  of  aluminum 
particle  size  is  used  to  ke /  the  fraction  burned  calculation, 

so  that  a  zero  value  indicates  that  this  option  will  not  be  used. 

•  Terms  factored  from  the  fraction  burned  and  slip  flow  equations 
are  calculated.  Factors  are  removed  from  the  calculation 
loops  to  reduce  execution  time. 

•  The  location  of  the  initial  head  end  of  the  motor  is  established. 
The  head  is  normally  the  physical  head  end,  but  is  changed  for 
special  cases  involving  igniter  location  and  fuse  configurations. 

•  The  following  groups  of  calculations  are  repeated  at  each  time- 
increment  of  the  problem.  The  igniter  flow,  if  any,  is  deter¬ 
mined  from  tabular  input,  an  outer-loop  convergence  factor  is 
calculated,  and  terms  factored  from  the  transient  deflection 
equations  are  prepared.  A  series  of  parameters  are  set  to 
initial  values. 
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•  The  particle  lag  terror  are  completed.  Since  these  values  are 
assumed  dependent  on  head  pressure,  they  are  constant  down 
the  port  and  are  not  recalculated  at  each  axial  station.  Gas 
properties  and  deflection  also  use  this  interpretation,  as 
follows . 

•  Input  values  of  thermochemistry  program  data  are  interpolated 
from  input  tables  as  a  function  of  the  fraction  ot  aluminum 
burned.  The  gas  properties  required  to  solve  the  equations  of 
motion,  such  as  the  ratio  of  specific  heats,  the  gas  constants, 
and  sonic  velocity,  are  calculated. 

•  The  transient  deflection  calculation  is  completed.  The  static 
shape  terms  from  tabular  input  were  placed  in  an  array  prior 
to  the  beginning  of  looping. 

•  The  following  calculations  are  made  for  each  axial  increment  in 
the  motor.  Trial  values  for  the  initial  time  step  are  used  from 
the  previous  axial  station;  all  other  trial  values  result 
from  the  previous  iteration  of  the  station.  Special  case  logic 
for  the  head-end  node,  ignition  options,  or  initial  time,  are 

not  shown  in  the  figure  or  discussed  here.  (However,  they 
should  be  recognized  in  the  program  listing  as  I.  EQ.  L, 
IGN.GT.O,  or  TIME.  EQ.  0.  branches.) 


•  The  port  geometry  is  calculated  from  the  web  burned,  burning 
rate,  and  the  input  options  of  a  circular  port,  two-dimensional 
port,  or  star  fuse  design.  Symmetry  is  assumed  in  all  cases. 

•  Base  burn  rate  for  the  propellant  at  the  axial  station,  as  a 
function  of  local  pressure,  is  found  from  an  input  table.  The 
input  table  can  be  keyed  to  use  log-log  interpolation.  The  ero¬ 
sive  contribution  is  then  calculated. 

•  Trap  logic  is  executed  to  deteimine  the  mass  addition  during 
tailoff  of  the  motor.  Some  local  conditions  are  calculated 
using  the  energy,  state,  and  continuity  equations. 

•  The  deflected  port  area  is  determined,  traps  are  executed  to 
prevent  crossing  Mach  1,  and  the  momentum  equation  is  used 
to  solve  for  pressure.  The  set  of  equations  at  each  axial 
station  is  iterated  with  logic  for  low- speed  flow,  high-speed  but 
subsonic  flow,  and  supersonic  flow. 

•  Looping  logic  allows  10  (or  a  larger  input  number)  iterations 
and  then  accepts  the  result  whether  converged  or  not.  The 
error  is  saved  for  printout  should  the  outer  loop  converge  when 
the  inner  loops  are  still  in  error. 
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•  The  geometry- flow  relations  used  to  determine  the  aerodynamic 
throat  are  checked,  and  logic  to  relocate  the  throat  node  is 
executed  if  necessary.  Any  changes  are  printed. 

•  The  rode  number  of  the  axial  increment  is  checked  to  branch  to 
the  next  step,  as  indicated  in  the  figure.  Throat  and  subsonic 
exits  must  branch  to  check  exit  constraints,  and  a  supersonic 
exit  bypasses  to  the  output  section.  All  other  locations  go  back 
to  calculate  the  next  axial  station. 

•  If  a  plenum  and  nozzle  are  modeled,  the  plenum  pressure  is 
calculated. 


An  updated  head  pressure  is  calculated  from  the  error  in  Mach 
number  or  a  term  created  to  indicate  an  overchoked  prior  cal¬ 
culation.  If  the  exit  pressure  is  less  than  ambient,  the  motor 
is  unchoked  and  head  pressure  is  iterated  to  force  p  ..  = 
^ambient’  ^  note  is  printed  as  the  unchoke  occurs.  (If  exit 

flow  is  supersonic,  the  exit  pressure  may  fall  below  ambient. 

It  is  assumed  that  separation  in  the  nozzle  does  not  occur.  ) 


•  If  the  throat  or  exit  constraint  is  not  satisfied,  loops  are 

counted  to  10  (or  a  larger  input  number).  From  10  to  20  loops, 
the  error  is  printed;  and  after  20  loops,  the  result  is  accepted 
if  the  error  in  head  pressure  is  less  than  5  percent  (normal 
runs  would  converge  exit  conditions  to  0.5  percent  in  Mach 
number).  Otherwise,  the  run  is  terminated.  The  printed  error 
messages  indicate  whether  the  calculations  are  oscillating  or 
are  overdamped  and  allow  the  user  to  determine  whether  the 
run  is  acceptable. 


•  A  check  is  made  on  unchoked  motors  to  determine  if  choking 
has  occured,  in  which  case  the  throat  is  located  and  motor 
conditions  are  reiterated  with  "choked"  logic. 


•  Parameters  along  the  motor  are  printed,  if  requested. 


•  Values  used  in  transient  terms  or  integrations  are  stored  for 
use  in  the  next  time  step. 

•  The  instantaneous  divergence  angle  from  the  throat  to  the  exit, 
as  well  as  thrust,  are  calculated.  Summary  terms  are  stored/ 
accumulated  until  thrust  decays  to  10  percent  of  maximum, 
defining  the  action  time  of  the  motor:  summary  parameters  are 
printed. 

•  The  L*  instability  warning  criteria  is  checked  until  a  warning  is 
printed  or  the  head  web  burns  out.  The  message  appears  ahead 
of  the  normal  output  for  the  time  increment. 
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•  Normal  output  parameters  arc  printed,  and  input  tor  a 
trajectory  program  is  printed  punched  if  requested. 

•  If  the  run  is  the  baseline  for  parametric  studies,  selected 
parameters  are  retained.  After  the  baseline  run,  the  stored 
values  are  used  to  calculate  sensitivity  coefficients.  The 
parameters  selected  are  those  observable  in  experimental 
firings,  namely,  head  pressure,  thrust,  exit  burn  rate,  and 
burn  rate  at  an  input  x-ray  location. 

•  A  series  of  terms  are  saved  for  continuing  calculations,  the 
time  increment  for  the  next  step  is  determined,  and  the  head 
pressure  is  guessed.  Unless  the  run  is  completed,  the  calcu¬ 
lations  proceed  to  the  next  time  increment  at  the  head  of  the 
motor. 

•  Results  of  parametric  studies  are  sequenced  and  output  is 
printed  if  this  option  is  used. 

The  indexing  scheme  employed  in  the  program  is  shown  in  Table 
III- 1  as  an  aid  in  following  the  details  of  the  logic  of  special  cases.  For 


L  J  I  NSTAR  N 


example,  if  I  is  greater  than  NSTAR,  the  supersonic  flow  conditions  in  the 
exit  cone  are  being  calculated. 


Symbol 

L 

I 

J 

NSTAR 

N 

NP 

NPAR 


TABLE  III-  1 
INDEX  SCHEME 

_ Desc  ription _ 

H~ad-end  node,  L  =  1  except  for  ignition  transient  option 

Axial  node  location  of  calculation 

Prior  axial  node,  usually  J  =  I  -1,  J  >0 

Throat  node 

Exit  node 

Propellant  number,  1  or  2 

Location  of  the  parameter  being  varied  in  parametric  studies 
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2.  SUBROUTINE  STRUCTURE 

The  nozzleless  ballistics  model  is  essentially  a  single  main  pro¬ 
gram  with  three  utility-type  subroutines,  as  illustrated  in  Figure  III- 3.  The 
subroutines  have  the  following  functions: 

•  DATA  IN  is  a  random  read  routine  that  interprets  input  of 
selected  variables.  DATA  IN  operates  under  the  control  of 
TABIN,  which  requests  the  reading  of  cards  with  integers  or 
real  data. 

•  TABIN  is  a  routine  which  prepares  input  of  tabular  data.  The 
addressing  for  the  tables  is  arranged,  indexes  are  initialized, 
and  limits  to  table  size  established.  Input  of  integers  and  non- 
tabular  program  inputs  are  simply  passed  back  to  the  main 
program  for  interpretation. 

•  TAB  is  a  routine  to  interpolate  (or  look  up)  values  from  the 
tabular  inputs.  Linear  interpolation  or  log-log  interpolation 
is  used  and  no  extrapolation  is  made;  rather,  the  value  at  the 
boundary  of  the  table  is  returned.  The  tables  are  assumed 
of  the  form 

z  =  f  (x,  y) 

with  a  call  statement  of  the  form 
CALL  TAB  (X,  Y,  Z,  No) 

where  No  is  a  table  number  assigned  in  the  MAIN  program. 

The  size,  shape,  type  of  interpolation,  and  order  of  table 
loading  are  established  by  input  of  the  tables  and  thus  are  not 
transmitted  to  the  main  program. 

Internally,  the  logic  is  established  to  search  the  table  by 
stepping  from  the  prior  point  to  the  new  one,  consistent  with 
an  iterative  program  with  highly  repetitious  call  statements, 
in  contrast  to  computing  indexes. 

Tabular  input  may  be  degenerate  in  x,  y,  or  both,  and  missing 
tables  are  interpreted  as  z  =  0  with  one  and  only  one  error 
message  printed  as  a  warning. 

All  other  functions  and  subroutines  used  by  the  program  are 
believed  to  be  standard  library  routines. 

The  groups  of  program  output  are  shown  in  Figure  III  -  3  to  illustrate 
the  relative  position  in  the  program  where  output  it  printed.  A  more  explicit 
example  is  given  in  Section  V. 
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Figure  III- 3  Model  Structur 
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SECTION  IV 

INPUT  REQUIREMENTS 


The  general  requirement  for  input  information  needed  to  run  the 
program  is  discussed  first  to  show  the  type  of  data  that  must  be  available. 
Second,  the  scheme  used  to  interpret  the  input  is  described,  and  finally, 
the  specific  word-by-word  requirements  are  given.  The  rationale,  in 
contrast  to  rules,  is  considered  in  Section  V. 

1.  INFORMATION  REQUIREMENTS 

Motor  geometry  description  must  include  the  initial  port  diameter 
and  grain  outer  diameter  given  as  a  function  of  motor  length.  An  exit  cone, 
if  any,  is  part  of  the  port  diameter  profile. 

Propellant  thermochemistr  ,_calculation  outputs  for  the  desired 
range  of  aluminum  fraction  burned,  F,  from  so-called  T*  runs,  are 
expected.  The  input  is  given  as  parameters  from  the  ISp  printouts  in  the 
units  of  standard  equilibrium  thermochemistry  programs  (Section  II-4). 

The  range  of  F  should  not  be  less  than  will  be  used  in  the  run  because 
extrapolation  is  not  allowed:  The  limiting  values  of  the  given  data  will  be 
used,  and  no  warning  is  provided. 

The  igniter  mass  flow  as  a  function  of  time  and  the  upstream 
spreading  rate  are  needed  if  an  option  using  these  parameters  is  requested. 

Grain  deflection  is  calculated  from  a  shape  function  from  stress- 
analysis  program  calculations  and  coefficients  for  the  transient  equation; 
both  are  required  if  deflection  is  to  be  included  in  the  model.  The  im¬ 
portance  of  deflection  to  the  simulation  results  varies  with  motor  design 
and  propellant  mechanical  properties,  so  a  decision  to  neglect  deflection 
must  be  made  in  each  case.  If  case  deflection  is  significant  relative  to 
a  rigid  case,  this  effect  should  be  included  in  stress-analysis  calculations. 

Erosive  burning  coefficients  fcr  the  propellant  and  base  burn  rate 
from  strand  tests  are  required.  Again,  the  input  must  cover  the  range  of 
pressures  encountered  in  the  run,  usually  from  peak  pressure  down  to 
abient  or  even  below  ambient  if  supersonic  exit  cones  are  simulated.  If 
needed,  the  data  must  be  extrapolated  by  the  user. 

Miscellaneous  constants  and  selection  or  logic  to  control  the  pro¬ 
gram  are  needed  to  complete  a  problem  definition,  and  are  discussed  in 
subsection  IV-3. 

2.  INPUT  SCHEME 

A  technique  allowing  maximum  flexibility  of  problem  input  has  been 
adopted.  This  scheme,  however,  should  be  understood  because  it  is 
important  in  running  multiple  cases  in  a  single  computer  access. 
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The  first  card  in  a  case  definition  is  a  header  card  for  identification 
of  the  problem,  Figure  IV-  1.  The  header  (1)  is  read  from  the  main  program 
and  is  immediate!  '  printed  on  the  output.  A  1  is  required  in  the  first  column 
(c  c  1);  absence  of  the  1  results  in  a  normal  termination  of  the  run. 

The  remaining  irput  is  controlled  from  the  subroutine  TABIN,  but  is 
actually  read  in  DATA  IN,  in  the  sequence  indicated  in  Figure  IV- 1.  All  of 
the  data  tables  (2)  are  read  first,  with  a  limit  of  50  tables  and  3000i0  words 
of  storage,  both  about  three  times  more  than  normal  input  for  one  case.  The 
tabular  input  is  assumed  to  continue  until  a  completion  code  (loc  4  =  -1)  is 
given.  The  input  in  the  remaining  two  sets  (3;i  and  (4)  is  then  interpreted  in 
the  main  program.  Integer  words  (3)  are  logic  controls  for  the  calculation, 
and  data  words  (4)  include  coefficients,  constants,  and  some  problem  logic 
control  words,  to  be  described  in  detail  in  subsection  4.3. 

The  input  is  transferred  to  working  storage  in  the  main  program  (5) 
if  non-zero  (nz)  values  are  given.  As  a  result,  only  changes  in  the  problem 
are  needed  in  successive  cases.  Also,  only  input  pertinent  to  the  case  and 
options  to  be  used  is  required,  so  unneeded  input  can  be  completely  ignored. 

Since  zero  input  is  ignored,  no  "filling"  is  required,  but  values  in 
the  cards  do  no  harm  if  desired  for  column  alignment,  etc.  One  limitation 
in  *he  method  is  that  options  can  not  be  "turned  off"  by  inputting  a  zero  in  a 
later  case. 

In  batch  runs,  any  constants  or  tables  may  be  replaced  by  new  input, 
except  with  zero  integers  or  constants  as  discussed.  The  tabular  inputs  are 
added  into  the  lists  and  storage  area,  and  the  internal  logic  is  readdressed 
to  the  new  input.  This  approach  has  two  implications: 

(i)  The  new  table  is  independent  of  the  old  in  size,  type,  and 
interpolation  method. 

(ii)  New  tables  add  to  the  total  number  of  tables  (50  limit)  and 
storage  (3000  words),  so  the  cumulative  number  in  batch 
runs  can  exceed  these  limits. 

3.  DETAILS  OF  INPUT 

In  this  subsection,  the  mechanics  of  card  format  and  table  input  are 
considered,  then  the  parameters  and  requirements  are  specified. 

a.  Format 

Three  card  input  formats  are  used  by  the  program, 

II,  71H  Header,  one  card 

4  (12,  14,  6X),  12,  28A  Logic 
6  (12,  F 10),  12  Data 


IS 
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Figure  IV- 1  Input  Groupings 
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For  either  logic  (integer)  words  or  data  cards,  the  input  is  interpreted  as 
(loc  for  location) 


1-2 

3-12 

1  3-  >4 

loc, 

input, 

loc 

\  alue 

12, 

14,  6X  or  FI  0, 

12 

The  input  cards  are  read  in  sets,  and  a  set  is  assumed  to  continue  until  a 
completion  sign  is  encountered.  Two  completion  signs  are  recognized: 

(i)  A  blank  card 

(ii)  A  negative  location,  i.e.,  location  =  -  1  in  any 
of  the  12  fields  on  the  card. 

Note  that  the  completion  sign  for  tabular  input  is  different  in  meaning  and 
form  (b4bb  -1  in  1^.  14). 

Rules  governing  the  interpretation  of  a  card  set  include: 

•  Any  number  of  cards  may  form  a  set. 

•  Any  order  of  locations  is  acceptable. 

•  Zero  values  of  input  are  ignored. 

•  If  locations  are  repeated,  the  last  value  of  input  is 
used  (overwrites). 

•  If  a  location  is  blank  but  data  is  given,  a  sequenced 
list  is  assumed,  so  the  following  two  cards  are 
equivalent: 

1  10.  2  11.  3  12.  -1 

1  10.  11.  12.  -1 

The  sequencing  does  not  follow  from  card  to  card,  so  the 
first  location  on  the  next  card  must  be  given. 

b.  T abula r  data 

Data  tables  are  assumed  of  the  form  z  =  f  (x,  y).  Tables  are 
termed  rectangular  if  the  number  of  x  arguments  in  the  rows  is  the  same  for 
each  of  the  columns  of  y  arguments.  A  nonrectangula r  table,  then,  has  a 
varying  number  of  *  arguments  in  the  columns. 
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The  card  sets  to  input  a  table  are  given  in  Table  IV-  1  with  the 
location  of  the  words  in  the  sets.  The  logic  (integer)  card  set  is  always 
required  and  must  include  a  table  number  and  a  completion  sign.  The 
number  of  y  columns  (NY)  in  location  5  and  x  rows  (NX)  in  location  6  mu'! 
be  given  if  more  than  one.  Location  3  is  used  to  key  log-log  interpolation 
if  desired,  and  a  (4)  value  in  location  4  keys  a  non r ec ta ngula r  table. 

Usually,  one  set  of  y  values,  one  set  of  x  values,  and  NY  sets 
of  /.  values  compose  a  table.  If  only  one  y  column  is  specified,  the  set  of  y 
values  is  omitted.  If  a  nonrectangular  table  is  specified,  the  number  of  \ 
values  in  each  column  must  be  given,  starting  in  location  7,  and  an  x- 
argument  set  is  required  for  each  y  argument. 

c.  Assignment  of  Input  Tables 

The  numbers  of  input  data  tables  are  assigned  in  the  main  pro¬ 
gram  due  to  compiled-in  constants  in  the  call  statement: 

CALL  TAB  (X,  Y,  Z,  NO) 

The  cor  responder,  e  of  input  data  and  the  table  numbers,  NO,  is 
given  in  Table  IV-2  with  the  engineering  symbol,  FORTRAN  symbol,  and 
assumed  arguments  of  the  tables.  In  Table  IV-2,  data  that  is  designated  as 
"optional"  is  needed  only  if  the  related  option  is  specified.  Data  that  is 
designated  as  required  is  almost  always  needed  since  zero  values  of  the 
parameter  (7)  are  rarely  meaningful;  the  possible  exception  is  the  list  of 
times  of  detailed  output  where  no  input  results  in  printout  of  all  time  steps. 

Symbols  in  Table  IV-2  a*-_-  separately  defined  in  Appendix  A. 

d.  Integer  Inputs 

Up  to  12  integer  words  may  be  input  to  control  program  options, 
as  showi  in  Table  IV-3.  Only  two  are  always  required,  the  number  of  nodes 
(increments  +1)  and  the  erosion  model  number  (olus  loc  (4)  =  NUM  (4)  =  -1  to 
indicate  the  end  of  data  tables). 

e.  Data  Input  Words 

Up  to  40  data  (real)  words  may  be  input  to  complete  the  input 
constants  and  control  options  in  the  program,  as  listed  in  Table  IV-4.  Note 
that  one  data  set  is  used,  regardless  of  the  number  of  cards,  as  discussed 
in  subsections  IV-2  and  IV-3. 

In  Table  IV-4,  data  required  for  most  runs  is  indicated.  Data 
related  to  optional  calculations  are  required  if  the  option  is  requested.  For 
example,  if  X2  is  given  in  Z(6)  to  "key"  use  of  two  propellents,  then  7(31), 
7(33),  7(35),  7(37)  are  needed  to  define  the  second  propellant  and  an 
erosive  burning  correlation.  In  addition,  the  base  burn  rate  should  be 
input  in  table  1  and  the  correlation  type  in  NUM  (3). 
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TABLE  IV  -  1 

INPUT  OF  TABULAR  DATA 


Input  for  Form  7.  -  f  (X,  Y) 


Input  T ype 

Input 

Desi  ription 

Integer  card  set 

Required 

Control  constants 

Data  card  set 

Requi red  if  NY  >1 

Omit  if  NY  =  1 

Y  values 

Data  card  sets 

1  set  if  rectangular 

NY  sets  ii  non rectangula r 

X  values 

Data  card  sets 

NY  sets  required 

7.  values 

Integer  Cards  4  (12,  14,  6x), 

12,  28A 

Locat  ion 

Input 

Desc  riptio"* 

1 

Req  aired 

Table  number  (sec  Table  IV-2) 

3 

Optional 

Key  to  LOG- LOG  interpolation  (+) 

4 

Optional 

Table  type,  nonrectangular  =  (  +  ) 

(Note  that  (-)  is  k.  y  to  end  table  inputs) 

5 

Optional 

Number  of  Y  (columns)  values,  required  if  NY  >1 

6 

Optional 

Number  of  X  (rows)  values,  required  if  NX  >1 

7  to  25 

Optional 

Number  of  X  values  in  2nd - 19th  Y  columns, 

nonrectangular  tables 

Data  Card  Sets  6  (12,  F10),  12 

Location 

Input 

Desc  ription 

1 

Optional 

List  of  X,  Y,  or  7.  values  sequenced  from 
Location  1 

Any 

Requi  red 

Completion  sign 

NOT  ES: 


1.  A  table  can  degenerate  to  one  point  in  X,  Y,  or  both. 

2.  Interpolation  is  linear  with  no  extrapolation:  values  at 
boundary  of  table  are  returned. 

3.  A  missing  table  causes  7.-0  and  one  ionly  one)  error 
message. 
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TABLE  IV -2 

INPUT  TABLE  ASSIGNMENTS 


No. 

F ur  ction 

Input 

Symbol 

z. 

F orm  z 

X 

(x,  y) 

y 

1 

Propellant  burn 
rate 

Required 

r 

o 

R0 

P 

PNUM 

2 

Igniter  mass  flow 

Optional 

w. 

ign 

W  (L) 

t 

3 

Motor  geometry 

Required 

d. 

l 

d 

o 

D 

DMAX 

X 

1  =  d. 

l 

2  =  d 

o 

4 

Axial  increment 
length 

Required 

dx 

dx, 

nom 

DX 

x/L 

— 

5 

Time  increments 

Requi  red 

dt 

DELT 

6 

— 

6 

Time  of  detail 
output 

Required 

— 

T  LIST 

Counter 

— 

7 

Nominal  deflection 

Optional 

— 

DEFLO 

X 

— 

8 

Chamber 
tempe  rature 

Required 

T 

c 

TC 

F 

— 

9 

Specific  heat  of 
products 

Required 

c 

p»  p 

CPP 

F 

— 

10 

Specific  heat  of 
gas 

Required 

c 

p.  g 

CPG 

F 

13 

Moles  of  gas 

Required 

— 

MOLG 

F 

— 

14 

Moles  of  products 

Required 

— 

MOLP 

F 

15 

Molecular  wt  of 
gas 

Required 

Mg 

MWTG 

F 

— 

1  1 

Time  of 
sensitivity  calc 

Optional 

— 

T  PART 

Counter 

— 

12 

List  of  parameters 

(parametric 

sensitivity) 

Optional 

- - 

NPAR 

Counter 

— 

16 

F raction  of  head 
burning 

Optior.al 

— 

HEADW 

t 

— 

IV  -  7 
« 
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TABLE  IV -  3 
INTEGER  INPUTS 


Loc  ation 

Symbol 

Input 

Desc  ription 

NUM  (h 

N 

Requi red 

Number  of  nodes  (increments  +1);  use 
10,  19,  or  28  (limited  by  output  format! 

NUM  (2 

KODE 

Required 

Erosion  model  number,  first  or  only 
propellant 

NUM  (3) 

KODE 

Optional 

Second  erosion  model  number 

NUM  (4) 

Required  =  -  1 

Keys  this  list 

NUM  (5) 

KEY 

Optional 

Code  word  requesting  calculation  of 
parametric  variations 

NUM  (6) 

IGN 

Optional 

Key  (  +  )  for  ignition  at  At  option 

NUM  (7) 

NSTAR 

Optional 

Throat  node  number,  a  slight  aid  to 
program  execution  speed 

NUM  (8) 

NXRAY 

Optional 

Node  number  of  local  r  output  for 
parametric  studies.  Not  used  without 
KEY 

NUM  (9) 

L 

Optional 

Noie  number  for  ignition  if  not  at  head 
of  motor 

NUM  (10) 

LOOPM 

Optional 

Number  of  loops  run  without  error 
message  (10  assumed  if  not  greater) 

NUM  (12) 

LIST 

Optional 

Key  (  +  )  to  output  for  trajectory 

program  inputs 
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TABLE  IV -4 
DATA  INPUT  WORDS 


Location 

Symbol 

Units 

Input 

Desc  ription 

*0) 

- 

TSTOP 

sec 

R  equi red 

Computation  time  if 
desired  less  than  burn¬ 
out  time 

Z  (2) 

TRISE 

sec 

Optional 

Approximate  time  of  fast 
pressure  rise 

Z  (3) 

- 

ALUM 

— 

Required 

Aluminum  fraction  in 
propellant,  required  if 
not  zero 

Z  (4) 

P 

BETA 

— 

Required 

Coefficient,  shear 
stress  calculation 

Z  (5) 

MU 

lb/in.  -  sec 

Required 

Viscosity  of  combustion 
products 

Z  (6) 

- 

X  2 

in. 

Optional 

Length  at  beginning  of 
second  propellant 

Z  (7) 

L 

XE 

in. 

Required 

Motor  length 

Z  (8) 

XPLEN 

in. 

Optional 

Plenum  length 

Z  (9) 

DPLEN 

in. 

Optional 

Plenum  diameter 

z  (10) 

DTH 

in. 

Optional 

Throat  diameter 

Z  ()  i) 

PPLEN 

psi 

Optional 

Initial  plenum  pressure 
(if  nz,  Z(8)  -  7(10) 
are  then  REQ'D) 

z  (12) 

WIDE 

in. 

Optional 

Width  of  star  sheet,  key 
to  star  geometry 
calculations 

Z  (13) 

HIGH 

in. 

Optional 

Motor  width,  key  to 
two-dimensional 
geometry  calculations 

Z  (14) 

STAR 

- 

Optional 

Number  of  points  in 
star  design 

Z  (15) 

VFZ 

in. /sec 

Optional 

Velocity  of  upstream 
ignition  propagation 
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TABLE  IV- 4  (Continued) 


Location 

Symbol 

Units 

Input 

Desc  r  iption 

Z  (16) 

WFZ 

lb/ sec 

Optional 

Mass  flow  due  to  fuse 
combustion  ( 7  (  1  S)  is 
then  REQ'D) 

Z  (17) 

EPSP 

Optional 

Term  to  determine 
port-to-plenum  pres¬ 
sure  rise  in  nozzled 
motor 

Z  (18) 

STEP 

- 

Optional 

Code  word,  restart 
parametric  variations 

Z  (19) 

DOME 

Optional 

If  negative,  key  to 
hemispherical  head 
geometry 

Z  (20) 

CONV 

- 

Optional 

Convergence  ratio, 
suggest  0.0005  (0.05% 
local,  0.1%  mid,  0.5% 
throat  accuracy) 

Z  (21; 

CHOKE 

- 

Optional 

Initial  conditions 
assumed  choked  if 
positive 

Z  (22) 

PSTART 

psia 

Required 

Initial  head  pressure  gues 

Z  (23) 

RHOAMB 

lb/in.3 

Optional 

Initial  gas  density  in 
port 

Z  (24) 

PAMB 

psia 

Optional 

Ambient  pressure 

Z  (2  5) 

FACTOR 

- 

Required 

Loop  gain  factor, 
suggest  1.! 

Z  (26) 

KD 

lb/in.3 

Optional 

Spring  constant  in 
transient  deflection 
calculation  ( 1  0k  -  50k) 

Z  (27) 

WN 

rad/ sec 

Optional 

Undamped  natural  fre¬ 
quency  in  transient 
deflection  and  key  to 
calculation  (40  -  100) 
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TABLE  IV -4  (Continued) 


Location 

Symbol 

Units 

Z  (28) 

PREF 

psia 

*  (29) 

ALPHAN 

- 

Z  (30) 

G02  or  BH 

- 

Z  (31) 

G02  or  BH 

- 

Z  (32), 

Z  (33) 

XKG  or  PR 

- 

Z  (34), 

Z  (35) 

XM 

- 

Z  (36), 

Z  (37) 

RHO 

lb/ in. 3 

Z  (38) 

DALUM 

mic  rons 

Z  (39) 

SALOX 

- 

Z  (40) 

SALUM 

- 

Input 

Dose  ription 

Optional 

Reference  pressure 
used  in  nominal  grain 
deflection  input 
(Table  A-2) 

Required 

Relative  value  of  L* 
correlation  to  70°F 
value 

Required 

First  coefficient  in 
erosive  burning  model, 
first  propellant 

Optional 

First  coefficient  in 
erosive  burning  model, 
second  propellant 

Required, 

Optional 

Second  coefficient,  as 
7(30)  -  7(31) 

Required, 

Optional 

Exponent,  as  7(30)  - 
7(31) 

Required, 

Optional 

Propellant  density 

Optional 

Aluminum  mean  particle 
diameter,  key  to  frac¬ 
tion  burned  calculations 

Optional 

Slip  ratio;  AlzO} 

Optional 

Slip  ratio,  aluminum 
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Similar  groupings  apply  for  fraction  burned  calculations 
Z(  iH)  -  Z(40)  and  for  a  motor  with  a  nozzle  Z(8)  -  Z f  1  1 ) ,  Z(17). 

A  special  comment  regarding  the  loop  gain  FACTOR  in  location 
Z  ( 2  5 )  is  required.  The  Mach  number  convergence  criteria  at  the  throat  has 
a  tendency  to  change  for  different  problem  definitions,  so  the  solution  can 
diverge  or  become  overdamped  and  require  excessive  iteration.  Output  is 
provided  to  indicate  the  type  of  difficulty  so  that  adjustment  can  be  made  in 
later  runs:  decreasing  the  input  FACTOR  damps  the  calculation. 

Parametric  studies  are  performed  by  perturbing  individual 
parameters  and  repeating  the  case  under  program  control.  From  differ¬ 
ences  between  results  at  specified  times,  a  sensitivity  parameter  is  calcu¬ 
lated  and  output;  the  output  is  necessarily  a  finite  difference  table  (in 
contrast  to  slopes).  Input  requirements  are: 

(1)  STEP  =  1.  to  key  logic 

(2)  TPART  in  table  11,  list  of  times  desired  in  sequence 

(3)  NPAR  in  table  12,  list  of  locations,  in  the  input  array  Z, 
of  the  desired  parameters  to  be  perturbed. 

The  number  of  output  parameters  and  times  are  limited  to  5  and  15, 
respectively,  by  dimensioning  of  array  rl(9,  15)  and  can  be  easily  modified. 
The  locations  must  be  input  as  real  variables  to  conform  to  the  table  input 
format. 
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SECTION  V 
SAMPLE  CASES 


The  sample  cases  given  in  this  section  are  intended  to  illustrate 
model  operation;  thus  they  are  not  specific  motor  designs. 

The  output  listing  of  a  sample  case  is  shown  in  Table  V-l  and  serves 
as  an  input  example  as  well,  since  all  input  is  output  on  the  listing.  The 
format  is  expanded  in  printing  but  the  output  represents  the  appearance  of  the 
input  cards.  The  words  LOGIC  ?nd  DATA  are  added  at  the  left  to  show  if 
integer  or  real  data  was  expected  by  the  random-read  subroutine  (for  diag¬ 
nosis  of  input  errors).  Comments  after  LOGIC  are  input  on  the  cards  and 
are  not  used  by  the  program. 

As  an  example  of  the  input  of  tabular  data,  the  third  table  from  the 
top  is  input  as  loc  =  1,  integer  =  3  to  identify  the  motor  geometry  table. 
Location  5  specifies  that  2  values  of  y  are  to  be  used,  and  location  6, 
because  it  follows  5  in  sequence/1' specifies  3  values  of  x.  A  completion  code 
(-1)  is  in  the  fourth  location  field.  The  two  values  of  y  are  sequenced  from 
location  1  in  the  card,  and  the  set  is  terminated  by  a  completion  code.  The 
three  x  values  form  a  set,  and  two  sets  of  z  finish  the  table.  The  data  in 
this  table  establish  a  port  diameter  (y  =  1)  tapering  from  0.621  inch  at  the 
head  to  0.64  inch  at  a  motor  length  of  39.24  inches,  then  flaring  to  1.46 
inches  at  a  length  of  40  inches.  The  grain  diameter  (y  =  2)  is  1.96  inches. 

The  last  two  groups  of  input  are  indicated  by  the  4  - 1  in  the 

third  field  at  the  end  of  the  table  input.  This  case  used  19  nodes  and  burn 
rate  model  1.  The  last  six  lines  of  DATA  are  the  remaining  input,  as 
one  set. 


The  comments  on  the  next  half-page  of  the  listing  summarize  the 
problem  setup  before  calculation  begins.  The  first  comment,  indicating  that 
no  table  11  was  input,  is  a  warning,  but  the  table  is  not  required  in  this  case. 
The  next  line  summarizes  the  erosion  model  for  propellant  1;  since  there  is 
only  one,  a  second  line  is  not  printed.  A  rough  summary  of  geometry  is 
printed  but  the  throat  is  assumed  at  the  exit  since  no  node  location  was  input 
(in  NUM  (7)).  The  list  of  options  shows  that  igniter  flow,  fraction  burned, 
and  grain  deflection  are  keyed. 

The  heading  at  the  top  of  the  next  page  indicates  the  axial  node 
locations  and  is  a  good  check  on  the  related  inputs.  As  calculations  within 
the  motor  begin,  the  throat  is  recognized  at  39.2  inches. 


(<) 


The  output  shows  "zero"  rather  than  the  assumed  location.  This 
choice  is  made  to  aid  in  diagnosing  errors  in  input  cards  (the  card 
"image"  is  returned). 


(The  reverse  is  blank) 


AFRPL-TR-73-20 


TABLE  V-l 

SAMPLE  CASE  OUTPUT  LISTING 


V-  5 


c0 
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6 

6 

4 

-1 

DATA 

0.000000 

-0 

1 • 000000 

-0 

2.030000 

-0 

3.000000 

-1 

DATA 

.109000 

-3 

1 .000000 

-0 

2.300000 

-0 

ICO.OOOOOO 

-l 

LOGIC 

1 

1 

l 

5 

1 

-0 

12 

DATA 

1 0. 0000 jO 

-0 

100. 000000 

-0 

300.000000 

-0 

500.000000 

*0 

1000.000000 

-0 

1500.000000  *0 

OAT  A 

2000.000300 

-0 

2500.000000 

-0 

3000.000009 

-0 

3500.000000 

-0 

4000.000000 

-0 

5000.000000  -1 

OAT  A 

.051000 

-3 

. 190000 

-0 

.322000 

-0 

.410000 

-0 

.5*0000 

-0 

•  670000  ••0 

DATA 

.  750300 

-0 

.830300 

•0 

.rfiOOOO 

-0 

.875000 

-0 

.945000 

-0 

1.070000  -1 

LOGIC 

i 

5 

2 

-0 

3 

-1 

DATA 

1.000000 

-3 

2.000000 

-1 

DATA 

O.OOODOU 

-3 

T9. 240000 

-0 

40.000000 

1.4600M 

-l 

DATA 

.621000 

-0 

.640000 

-0 

-l 

JAT  A 

1.960000 

-0 

1.960000 

-0 

1.960000 

-1 

LOGIC 

5 

6 

3 

-1 

-Mb 

A 

data 

DATA 

LOGIC 

0. oooooo 
.010-900 

-0 

-3 

.100000 

.050000 

-0 

-0 

•  ioo£> 

.100003 

-1 

-l 

rSepr 
1  best 

Irom 

0,dva‘ubUi£Et. 

9 

) 

2 

6 

2 

-l 

A-  ■*" 

DATA 

0.000000 

-3 

.030000 

-l 

DATA 

. 550000 

-0 

LOGIC 

4 

5 

6 

-1 

OAT  A 

. 190000 

-c 

.400000 

-6 

.890001 
. 90006# 

-0 

•900000 

-0 

.940000 

-0 

.941000  M 

DATA 

1.900000 

-3 

.900000 

-0 

-0 

.555000 

-0 

.355000 

-0 

.344000  -1 

LOOK 

7 

6 

1 

-1 

DATA 

11 .390000 

•0 

19.240000 

-0 

40.000000 

-1 

DATA 

0. 000003 

-0 

.010000 

-0 

LOGIC 

ibr  lp: 

8 

6 

5 

-1 

DATA 

0.000300 

-3 

.250000 

-0 

.530000 

-0 

.750000 

-0 

1.000000 

-I 

DATA 

22  70. OuODOO 

-0 

2oS0. 003000 

-0 

2966.000000 

-0 

3214.000000 

-0 

3377.000000 

LOGIC 

0 

A 

5 

-1 

DATA 

U.  0000.10 

-0 

.250000 

-0 

.500003 

-0 

•TfOOOO 

-0 

1.000000 

DATA 

9.  782100 

-3 

10.417500 

-0 

10.965900 

-0 

11.426300 

-0 

U.  790000 

"1 

LOGIC 

10 

6 

5 

-1 

DATA 

0. 000000 

-3 

.250000 

-0 

-9 

.750000 

-0 

1.00  OOM 
9.119100 

DATA 

10.255703 

-0 

10.060400 

-0 

-0 

9.457700 

-0 

*1 

LOGIC 

13 

6 

5 

-l 

DATA 

D. 000 J  00 

-0 

•250000 

-0 

.500000 

-0 

.750000 

-0 

1.000000 

0«TA 

3.424000 

-3 

J. 434000 

-3 

1*440906 

-0 

S.5O5O0Q 

-0 

6.5*7000 

LOGIC 

14 

6 

5 

-1 

DAT  A 

0.000300 

-0 

.  2  50  GOO 

-0 

.500003 

-0 

.750000 

-0 

1.000000 

OAT  A 

4.165099 

-0 

4.062000 

-0 

4.014000 

*0 

3.957000 

-0 

3.19*000 

“1 

LOGIC 

15 

6 

5 

-1 

DATA 

0. JOOOOO 

-0 

.250000 

-0 

.500000 

-0 

.750000 

-0 

1 . OOOOOO 

DATA 

23.360300 

o 

22.010000 

-0 

20.6200M 

t 

-0 

19.35000 0 

-0 

16.420000 

LOGIC 

19 

•0 

l 

4 

-l 

-1 

DATA 

10. 030300 

4 

53.000000 

-0 

.030004 

7 

40.000000 

3 

.200000 

-0 

DATA 

.0005  00 

-3 

1 .OOOOUC 

-0 

1100.000003 

-0 

.000100 

-0 

1 4. 900000 

-0 

.600000  -0 

DATA 

’6 

6300. OOOOOC 

-0 

40.000000 

-0 

1000.000000 

-0 

1.000000 

-0 

DATA 

It) 

20.0000  DO 

-0 

.947000 

-0 

.554300 

-0 

DATA 

H 

l .000000 

27 

.689000 

-0 

DATA 

ID 

.628000 

1? 

.  1  78400 

36 

.064000 

-1 

REFER  f  N 

:e  fit 

UN CORRELATED 

TABLE  NUMBER  11  0 

EROSION 

NODE! 

N  l  1  EOR  PRUPt  L  L  ANT 

1  REO 

- 

0440  K(  - 

.6280 

Kl  *  .1704 

EXP  •  1.00 

TO  LENGTH  ■  0.0  IN. 

NOMINAL 

r.r  :he 

THY 

HF  AD 

T  HR U AT 

EXIT 

LENGTH 

43.03 

40.00 

|T<NER  III  AN 

.62 

1  46 

1.46 

OUT  P R  Oj  AM 

1.96 

1  .9<> 

1.96 

OETAILEO  UFSIGN  INPUT  IN  TABLE  3  10  U$F  1?  NODES  SPACEO  AS  INPUT  IN  T  ABU  * 


CQNTR  Ql  NURDS  IN  INPjT  IN  1N-IFRG 
P  A  9  AH  *  MIC  VARIATION  AFY 
I.NITKN  INPUT  OT  OPT  I JN 
2-  *  I M '  NS I  'NAL  ut-  AIN  Kc* 
N-EiUnT  STAR  GRAIN  KEY 
fuse  I  ON  It  ION  CENT  ROC 
HE  AO- E  NO  'TflNf  DlSIGN 
PLf  NUH“TH^UAT  <  AlC  •  P  • 

►  PACT  1  ON  MU-Ni O  K\ V  0  » 
G-*  A T N  IJE  F  L  E  L  T  | oh  *E  Y  *N 


S  *£V| 

0  XRAY  NUDE  AC  IS 

0 

Q.OUO  IHIGHI 
0.030  (MlOEl 
0.0  I  RAT  El 
I.G 


20.0  MICRONS 
AO  .0 


Y  -  4 
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p.-ssu*£  thrust  Finn 

T|NF  HF  A  )  Ult  PL6N 

p*  rr*  ISP  imp  »jm  m  OIFl 

0.00 

8.00 

Lb. 00 

MEb 

20.00 

burned 

24.00 

/  l 6 NOTH 

28.00  32. uu 

36.00 

38.47 

40.00 

1008 

f  HR  DAT  t  tiC  4  r  tu«|  AT  NflOl  10  AT  LENGTH  39.2  IN,  SlCPE  - 

1.34 

.01  liei.'i  4?. 4  14.0  4)5.  2.U90 

553.5  Iht.b  20(3.1  2.  .010  .  000 

10  -2.88UI -02  2.44166*0)  3.25796-01  1.00006*00 

11  -l.74)6'-02  2.44146*03  3.25796-31  1.00006*30 

•  003 

.003 

.004 

.004 

.004 

•  004 

.004 

.005 

.005 

.001 

5 

.02  2441.)  TJ.4  14.0  1041.  4.652 

1129.9  1979.6  214.6  li.  .15?  .001 

.013 

.014 

.015 

.016 

.018 

.019 

.021 

•  023 

.025 

.006 

12 

.34  2613.1  1  ?  5 . 3  14.0  1 40b.  6.545 

1131.5  2J40.3  214.  d  35.  .171  .009 

.029 

.030 

.035 

.030 

.041 

.046 

•  061 

.057 

•  Ot; 

•  017 

1 

.07  2226.  7  143.9  14.0  1508.  7.096 

957.0  1702.5  2  12.5  7b.  .158  .016 

.060 

.053 

.062 

•  067 

.074 

.082 

.092 

,  103 

.113 

.03  7 

5 

1*0 

G 

2 

T- 

A 

M 

R 

RHO 

l 

0. 

0. 

On 

!■  71!  126*03 

4.76146-01 

0. 

7. 04706-01 

4.69146-04 

2 

0. 

o. 3ttd?t- Jl 

*ilOU2E*03 

1.72606*03 

4.80096-01 

4.55066-01 

7«  06356-01 

4.67906-04 

3 

T.  04  8  76- 03 

1.267)6* 33 

4.21 65F*01 

1.  71  106*03 

4.92)36-01 

9.25056-01 

7.4291E-U1 

4.64126-04 

4 

1.06696-02 

1.89526*03 

6.  3  9616  *03 

1.  Ob49E  *03 

5.09036-01 

1.43166*00 

7.92626-01 

4. 5755t-04 

5 

2.  2054F-O2 

2.52336*00 

8.  70I3E  *03 

1.64586*03 

5.30006-01 

l.9847E*00 

8.54256-01 

4.47696-04 

6 

3.  30  66*  -  02 

2 . 6  3b6F  *90 

9, 9 248E  *0 3 

1.62076 *03 

5.42256-01 

2.28236*00 

4.  89036-01 

4.41346-04 

7 

3.32436-02 

3. i  493E*  30 

1.  12086  *04 

1. 59126  *C3 

5 • 5561 E -01 

2.59646*00 

9.2;  366-01 

4.33886-04 

A 

4.  7568F-J2 

3.45086*00 

l.  25696  *04 

1. 5 568E  *C3 

5.702/E-01 

2.92846*00 

9. 66206-01 

4.25176-04 

9 

>.  795  7F*  02 

3,  7  7056*  30 

1.41296*04 

1.51686  *03 

5. S6i3E-01 

3.28036*00 

1.00856*00 

4.  15016-04 

10 

6. 9bJ;t-02 

4.0  779E  *00 

1.  5619F  *04 

1.47016*03 

6.03046-01 

3.6  5  3  76  *  00 

1.05286*00 

4.03I6E-04 

II 

0. 266?*- J2 

4.  3  81 3 E ♦ DU 

1.  7 1 76E  *04 

1.41556*03 

6.21076-01 

4.05056*00 

1.09196*00 

1.89456-04 

12 

9. 74)oF -u2 

4.677bE*00 

1.9  J67E*Q4 

1.35106*03 

6. 44476-01 

4.47316*00 

1.14766*00 

3.  73016-04 

13 

1.131 4F-JI 

4 . 94356*00 

2.  1  5  766  *04 

1.27606*03 

6.67206-01 

4.92  306*00 

1.19506*00 

3.53026-04 

14 

1.234  OF -  31 

6.140)6*00 

2.  3953F  *04 

1.1906c *03 

6. 097Qt-Ol 

5.39956*00 

1.23566*00 

3.31  786-04 

15 

I »  47  J2F • 01 

6.3345EOO 

2.  7046E  *04 

1.08396*03 

7.14996-01 

5.90396*00 

1.27796*00 

J.  04606-04 

16 

1.6I28E-J1 

5.43406*30 

2.96746*04 

9.98796*02 

7.32986-01 

6.22956*30 

1.30506*00 

2.02666-04 

17 

l.  (32  356-lU 

5.50956*00 

3.3  7  706*04 

8.78106*02 

7.54786-01 

6.56856*00 

1.33696*00 

2.51706-04 

Id 

>.  17276-01 

5.50716*00 

3.93516*04 

7.37066*0 2 

7. 75166-01 

6.78708404 

1.36066*00 

2.16106-04 

n 

9.2542fc-J2 

2.3  9466*00 

7.  1 557E  *04 

1.42471  *0? 

1.96406*00 

6.98  0  26400 

6.52946-01 

5.16536-05 

PR=SSU46 

THRUST 

H 

.On 

W66 

BURNED 

/  LENGTH 

T  IMF 

H*- 

A  ) 

F* 

1  T 

PL* 

•  N 

0.03 

0.90 

16.90 

20.00 

24.00 

28.00 

32.00 

36.00 

38.47 

40.00 

LOOP 

• 

P 

T  * 

ISP 

IMP 

S 

UM  M 

OEFL 

.11 

1711 

•  2 

14?, 

.5 

14. 

.3 

1450. 

6. 

900 

.079 

•  083 

.o§6 

.105 

.116 

.128 

.143 

.150 

.170 

.06  3 

0 

717 

.1 

1297 

203. 

.6 

m 

.628 

.018 

..1 

110  3 

.  5 

125, 

.5 

14, 

,0 

321. 

6. 

499 

.143 

.150 

•  172 

.187 

.204 

.223 

.244 

•  266 

.213 

.  122 

6 

Hb  f 

•* 

8’.? 

.6 

201, 

,2 

272 

1 

.302 

.011 

.  31 

7  7. 

.2 

104. 

.0 

14, 

.0 

139. 

5. 

758 

.198 

*203 

.233 

•  251 

•  271 

.294 

•  319 

.344 

.362 

.  165 

2 

329 

.4 

5  79 

.9 

197. 

,9 

395 

1 

.915 

.005 

-4| 

5H5 

.9 

92, 

.  3 

14, 

,0 

00  0. 

5. 

109 

.245 

.252 

.283 

.303 

.326 

.350 

.377 

.404 

.423 

.204 

4 

252 

.d 

445 

.  1 

194. 

,3 

503 

2 

.4  62 

.003 

.51 

47b 

.  J 

01. 

,4 

14. 

.0 

910. 

4. 

772 

.26  7 

.298 

•  348 

.372 

.■*■98 

•  426 

.453 

.472 

.235 

4 

235 

.  7 

>62 

.  1 

190. 

7 

599 

2 

.  760 

.002 

.6  l 

41?, 

.  4 

89, 

.2 

14. 

0 

054. 

4. 

659 

.325 

.332 

.365 

.388 

.413 

.440 

.468 

.495 

.513 

.252 

5 

1  74 

.  a 

307 

.0 

183. 

,4 

68  7 

3 

.432 

.001 

.71 

J5J 

.5 

7), 

.3 

1  4. 

0 

7  7  7. 

4. 

170 

.361 

•  367 

•  401 

.424 

.449 

.477 

.505 

.533 

.650 

.252 

5 

151 

.  4 

?  7  J 

.0 

13b. 

1 

760 

9 

.  874 

.001 

.8  1 

112, 

,  !, 

75. 

9 

1  4. 

0 

734. 

4. 

02  7 

.195 

.401 

.431 

.457 

•  483 

.519 

.539 

.  566 

.583 

.252 

5 

130 

.  1 

24? 

.2 

182 . 

4 

844, 

4 

.284 

.001 

.  9  l 

28), 

77. 

8 

14. 

0 

699. 

3. 

916 

.42  7 

•  433 

•  466 

.487 

•  514 

.541 

.570 

.597 

.612 

.252 

6 

125 

.6 

221 

.  1 

173. 

,6 

916, 

4. 

.601 

,000 

TAD 

V 

p 

A 

w 

R 

RHO 

1 

0. 

0. 

0. 

2.61116*02 

1.85346*00 

0. 

1.01236-01 

7.40906-05 

2 

0. 

1.38316-01 

2.012  6E  *01 

2.59716*02 

1.84766*00 

3.72  2  56  -01 

3.00466-01 

7.451  76-05 

3 

0. 

2.71966-01 

5.61366*03 

2.  55568  ♦  C2 

1 . 88936*00 

7.44436-01 

2.98146-01 

7.14146-05 

4 

0. 

1.94O6E-0L 

8.34006  *0) 

2.40906*02 

1.95476*00 

1.1  1  7  76*00 

2.94386-01 

/•  1  c406  -05 

5 

0. 

S.0269E-01 

1. 09936  *04 

2.40116*02 

2 • 0504 1  *09 

1.49346*00 

2. 89346-01 

6. 9291 t “05 

6 

0. 

5.5070E-01 

1.22856  *04 

2.  1 503  E «  02 

2. 10056*00 

1.68236*00 

?. 06-96-01 

6.  7V30E-05 

7 

J. 

5.  9473E-01 

1. 35616*04 

?.  29  55  F  *C2 

2.17286*00 

1.87226*00 

2.83166-01 

6.64  586-05 

8 

9. 

6. 35216-01 

1.40306*04 

2.2365  .*02 

2.24286*00 

2.06416*00 

2.00456-01 

6.43726-05 

9 

0. 

6. 72906-01 

1.  61466*04 

2.1  72c6*C? 

2.31  726*09 

2.25916*00 

2 .81 01 F  —01 

6.3 155* -05 

10 

9. 

7.0B44F-0I 

1. 75106*04 

2.  1  033  E  *0? 

2 . 3946F  *U0 

2.45  766*00 

2.81 2o6 -  01 

6.  Ue  Jfc-  05 

1  1 

J. 

7.4  IHjE-OI 

1 . 8976F  *04 

c.  J?  76  6  *C2 

2.47516*09 

2.66 016*00 

2. 61126-01 

5.  92iLt -05 

12 

0. 

7.  72926-01 

2.0551:  *94 

l.  94456*02 

2. 5599fc*0J 

2.06576*00 

2.8  1476-01 

5.u94l 6-05 

1  ) 

0. 

o. 02796-01 

2.2  J)2F  *04 

1.8511 F *02 

2.64306*00 

3.07426  *00 

2. 75256-01 

5.44726-05 

14 

0. 

H. 32876-01 

2.4496E  *04 

l.  74166  *02 

2.72191*00 

3.20506*00 

2. 77326-01 

5.1 5206-05 

15 

0. 

».fciboe-oi 

2. 72346*04 

1 . 60 B5E  ♦  C? 

?. 00006*00 

3.49716*00 

2.7396F-01 

4.79 J9L-05 

16 

). 

A.  78776-01 

2.94  79- * 04 

1.50546*02 

?.04d6f  *00 

3.62  796*00 

2. 70616-01 

4.51736-05 

1  7 

J. 

d. 9  59riF -01 

3.28036  *04 

1.  3  b  4?  E  *  C2 

2.(9306*00 

3.75  766*00 

2.64476-01 

4.  l40ik -05 

li 

9. C  95  76-01 

3.  d l 78F  *u4 

1. 16586*02 

2.91906*00 

3.03656*00 

2. 549  36-01 

3.61416-05 

19 

0. 

P. 77546-01 

4. 72556  *64 

d.  68576  *C1 

3.0  726*00 

3  .  0  3  6  56*00 

2.24756-01 

2.8127L-05 

V  -  5 

•> 

*  . 
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TmAuST  F  (  On  ttUPNfl)  I  U«ofH 


TIME 

HI  1 

i  > 

f  3  1  t 

Pll 

0.0) 

1.00 

16.30 

N. 

c 

O  1 

24.00 

26.00 

12.00 

It  .00 

16.4/ 

40.00 

LOOP 

P* 

IS*’ 

IN.* 

SUN  m 

Df  fi 

1.01 

261  . 

1 

Bb, 

,  4 

14. 

0 

6M. 

1. 

•  37 

•  456 

.46) 

.494 

.516 

•  54  J 

.9/0 

.999 

.625 

.640 

•  292 

6 

1  w 

.6 

2  )4. 

a 

17*. 

ft 

964. 

5. 

.  06  4 

.00) 

1.11 

24  1 . 

/ 

10. 

.  b 

l  4  ■ 

0 

64  5. 

3. 

,  15/ 

.417 

.492 

.522 

.5  44 

.9/0 

.996 

.  626 

.692 

.640 

.292 

6 

1  13 

.  1 

r#%. 

1 

1  71. 

A 

1050. 

5. 

,  449 

.003 

14  2. 

45*71 

-  33 

1.2  1 

21,. 

& 

90. 

.  1 

14. 

a 

47*. 

). 

150 

.s:  6 

.520 

.549 

.571 

.596 

.62* 

.691 

.661 

.660 

.252 

6 

<9 

.  H 

l  12. 

t 

1  73. 

6 

mi. 

5  , 

,  #04 

0, 

,003 

1.  il 

17b. 

J 

60. 

.  e> 

14. 

0 

462. 

2. 

76« 

.542 

.946 

•  $?4 

.595 

.619 

•  649 

•  66.' 

.661 

.660 

.  292 

1 

tl  J 

.  6 

141  . 

'i 

165  . 

V 

1163. 

6. 

Ill 

0 

.JQj 

1.41 

1 5  J. 

A 

TO. 

.0 

14. 

0 

395. 

2. 

407 

.964 

.570 

.597 

.617 

.640 

.66) 

.662 

.661 

.640 

.292 

9 

/O 

•  3 

l?». 

1 

164. 

1 

1  236. 

6. 

,  369 

0, 

.000 

10  -  I  .9/*‘fF  -  \ 
M  -l.O^Af'T 


L«  tfAHNlNG  AT  17.28** 


si  m. 

4 

57.6  14.0 

117.  1, 

.96) 

’OMf 

57.6 

►  -1  .SVPjF 

101. 

-01 

5  159.0  1241, 

6.569 

0.000 

61  Ml. 

1 

-7.  '  14.0 

2S5.  1. 

645 

.606 

4  7. 0 

>4. 

1  1  4.0  1270. 

*.770 

0.000 

71  61. 

/, 

39.0  14.0 

199.  1. 

343 

•  626 

39. 

6b. 

7  h».)  1292. 

6.920 

0.000 

*** 

61? 

.616 

.67-6 

.66) 

.662 

.661 

.640 

.292 

♦ 

j-*  . 

.611 

.6)6 

.65) 

.664 

.66) 

•  662 

.661 

.660 

.252 

U(* 

•  691 

.669 

.664 

.66) 

•  662 

.661 

.660 

•  252 

b  -l.664bt-'J) 

ArffVAf.f  l*UOS»  A^D  ISP  AVf  //3. 
1AM-J1  ThAOST  IS  1506.  LB. 


1.61  5,.  2  ?0.6  U.O 

?b.b  V).  ?  1  35.  7  1  309. 
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The  normal  output  at  each  time  increment  is  printed  in  two  lines 
under  the  heading.  The  error  messages  headed  10  and  1  1  indicate  that  the 
program  was  slow  to  converge  the  calculation  at  0.02  second,  and  was  over- 
damped,  Detail  output  was  requested  at  0.1  second  in  table  6,  the  first 
table  loaded,  and  printed  before  the  0.1  l -second  normal  output 

Calculation  continues  to  1.51  seconds.  The  L’:'  warning  indicates 
potential  instability  at  that  time.  Only  one  message  is  printed,  but  the 
nozzleless  motor  moves  toward  the  unstable  regior.  for  the  completion  of 
the  run. 


Between  1.41  and  1.5  seconds,  two  error  messages  are  printed. 
In  this  case,  an  error  of  0.2  and  0.1+  percent  in  velocity  was  accepted, 
where  0.1  percent  was  requested,  at  nodes  10  and  11,  respectively. 

Motor  action  time  is  sensed  at  1.81  seconds  as  the  thrust  drops  to 
.10  percent  of  peak  thrust.  A  summary  printout  of  average  thrust  and  IBp 
is  printed  ahead  of  the  normal  output  for  that  time  increment.  The  throat 
unchokes  prior  to  the  2.Cl-second  increment,  and  burnout  is  sensed  after 
2.11  seconds  as  head  pressure  falls  to  within  2  psi  of  ambient  pressure. 

The  next  three  sample  cases  (Table  V-2)  were  run  as  a  single 
machine  access  to  illustrate  the  input  of  changes  after  one  case  is  run. 
Also,  the  input  data  is  simplified  to  the  parameters  required  for  the 
specific  examples. 

The  first  case  in  the  batch  is  a  5-  by  40-inch  motor  with  a  1-inch 
bore  and  a  cone  to  2  inches  in  the  last  2  inches  of  the  motor  as  input  in 
table  3.  The  ignition  time  option  is  keyed  by  NUM(6)  =  1  input  in  the  card 
labeled  CONTROLS.  The  ignition  time  is  input  as  the  first  increment  in 
the  _\t  table,  table  5.  An  igniter  mass  flow  is  not  normally  consistent 
with  starting  with  the  motor  "on"  at  equilibrium,  so  is  not  input;  the  note 
of  the  missing  table  is  printed  on  the  second  page. 

The  summary  of  control  words  verifies  that  the  ignition  time 
option  is  set  and  that  the  fraction  of  aluminum  burned  calculations  are  set. 
Note  that  grain  deflection  is  not  to  be  considered,  and  also  that  no  data 
was  given  (Z(26)  -  Z(28)  or  table  7). 

At  this  point,  it  is  convenient  to  note  that  the  user  may  utilize  the 
"UNCORRELATED  TABLE"  notes  that  signal  a  missing  table  (11,  7,  2, 
and  16  in  this  example)  or  may  input  a  1  x  1  table  with  Z  =  0  to  avoid  the 
notes.  The  resulting  calculation  is  unchanged.  As  calculation  begins, 
the  remaining  two  missing  tables  are  inspected  and  the  throat  is  located 
at  node  1  7. 

The  first  printout  time  occurs  at  0.03  second,  as  expected.  The 
calculation  assumes  no  capacitance  term  in  the  continuity  equation 
(  9  PA  0)  with  this  option,  but  requires  input  knowledge  of  the  time  to 
9  t  peak  pressure. 


(  I  he  reverse  is  blank  ) 
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1.491 

1.717 

1.902 

1  .960 

2.013 

2.060 

2.098 

2.118 

2.114 

1.983 

10 

J.tO  172.2  M9.1  13.0  21*3.  13.002 

43.1  13*. 4  1 6M • 7  9*0  3.  33.03*  0.0 

!.*«0 

l.*2* 

L.  *4  3 

2.001 

2.093 

2.099 

2.13* 

2.194 

2.149 

2.014 

10 

3.70  144.2  13.4  14.0  2204.  13.103 

M3. 3  144.3  14M.3  4023.  34.343  0.0 

1.733 

Will 

1.  *4 

2.041 

2.093 

2.13* 

2.W4 

2.190 

2.1*4 

2.044 

10 

3.40  149.3  *1.7  14.0  2211.  13.201 

Ml. 7  193.3  144.0  4294.  33.4*0  0.0 

1.779 

1*912 

2.029 

2.0*1 

2.132 

2.176 

2.2  M 

2.226 

2.219 

2.074 

10 

3.90  160. b  MO.O  14.0  2229.  13.2** 

40. 0  190.3  167.6  694*.  37.009  0.0 

1  .  *1 4 

1.991 

2.0*9 

2.120 

2.171 

2.214 

2.24* 

2.2*1 

2.293 

2.»Cl 

10 

9.U0  197.2  7*. 9  14.0  2230.  13.3*2 

71.9  137.*  147.2  4490.  38.339  0.0 

1.894 

1.9*4 

2.109 

2.199 

2.209 

2.291 

2.2*4 

2.296 

2.207 

2.13? 

10 

9.10  153.8  77.1  19.0  2291.  1  3.4*9 

77*1  1  39.2  146.9  4914.  39.6*3  0.0 

1.8*7 

2.034 

2.144 

2.19* 

2.247 

2.288 

2.320 

2.331 

2.321 

2.160 

10 

4.20  190.7  79.4  14.0  2261.  13.976 

79.6  1  32.7  164.9  7140.  91  .036  0.0 

1.937 

2.073 

2.1*3 

2.21* 

2.2*4 

2.329 

2.396 

2.369 

2.394 

2.188 

10 

4.30  147.6  79.3  14.0  2271.  13.644 

79.3  1  30.3  1 46 .2  7366.  42.398  0.0 

1  .*7  4 

2.113 

2*221 

2.274 

2.321 

2.342 

2.391 

2.399 

2.3*7 

2.216 

10 

4.40  149.7  73.0  14.0  2280.  13.74* 

73. 0  128.0  169.8  7994.  43.76*  0.0 

2.019 

2.191 

2.29* 

2.111 

2.198 

2.398 

2*426 

2.433 

2.420 

2.244 

10 

4.90  119.0  97.7  14.0  1792.  11.147 

2.092 

2.1** 

2.2*4 

2.146 

2.392 

2.431 

2.43* 

2.438 

2.438 

.2.270 

11 

.7 

101 

.3 

160.8  7798. 

49.013 

0 

0 

I*U 

G 

V 

f 

A 

N 

R 

RMO 

1 

0.0 

0.0 

0.0 

*•17*26 

01 

1.446*6 

01 

0.0 

2.90446-01 

2.408  76-05 

2 

0.0 

5.10376-02 

3.1*726 

03 

*•  1 1 90E 

01 

1.94196 

01 

1.10996 

00 

2.89676-01 

2.39316-05 

3 

0.0 

9.74781-02 

6*20096 

03 

7.99*16 

01 

1.63936 

01 

2.24796 

00 

2. *6926-01 

2.34936-05 

4 

0.0 

1  .40091-01 

9.17936 

03 

7.70991 

01 

1.72906 

01 

3.40*06 

00 

2. *2336-01 

2  28006-05 

9 

0.0 

1  .01146-01 

i. 24026 

04 

7.39046 

01 

1. 79096 

01 

4.972*6 

00 

2.76236-01 

2  1*276-05 

6 

0.0 

1 • 99836-01 

1.40426 

04 

7.13436 

01 

1.629*6 

01 

9.19996 

00 

2.72336-01 

* . 12406-05 

7 

0.0 

2. 17876-01 

1.9*246 

04 

6.  *9326 

01 

1.86*16 

01 

5.73*16 

00 

2.67886-01 

2.05736-05 

8 

0.0 

2 . 39  336-01 

1.77936 

04 

*•*1936 

01 

1.90426 

01 

*.31786 

00 

2,62666-01 

1.98106-09 

9 

0.0 

2 .32  316-01 

1.99216 

04 

4.293*6 

01 

1*997*6 

01 

6.89346 

00 

2.56486-01 

1.89286-05 

10 

0.0 

2.493*6-01 

2.292  76 

04 

9.90746 

01 

1.9*396 

01 

7.49726 

00 

2.40*36-01 

1.78716-05 

11 

0.0 

2. *0916-01 

2. 497U 

04 

9*41*96 

01 

1.96396 

01 

7.77646 

00 

2.42456-01 

1.70856-05 

12 

0.0 

2.12106-01 

2.47916 

04 

9.9*776 

01 

1.9*396 

01 

7.80936 

30 

2.42306-01 

1.70006-05 

1  i 

0.0 

2.82226-01 

2.4*156 

04 

9.30406 

01 

1.9619E 

01 

7.8L2B6 

00 

2.422  76-01 

1.69966-05 

14 

1  .0 

2. *2 346-01 

2.4*326 

04 

5.58446 

01 

1.96336 

01 

7.81606 

00 

2.42246-01 

1  .69926-05 

13 

U.v' 

2  .82456-01 

2 .4*4*6 

04 

3.58296 

01 

1.96356 

01 

7.81916 

00 

2.4.206-01 

1.69886 -05 

lb 

u.u 

2. *2516-01 

2.48566 

04 

3.58216 

01 

1.96136 

01 

7.82076 

00 

2.42196-01 

1 .69066-05 

1 7 

0.0 

2. *2566-01 

2.40636 

04 

9.50146 

01 

1.96356 

01 

7.82226 

00 

2.42176-01 

1 .69046-09 

18 

0.0 

2.82906-01 

2.49206 

04 

5.57236 

01 

1.96356 

01 

7. 83 1 66 

00 

2.41996-01 

1.69606-05 

14 

u.o 

3.27996-01 

3.76676 

04 

4.09196 

01 

1.74226 

01 

*•05676 

00 

2,08856-01 

1.30116-03 
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TM«U$T  ►LOW 


WEB  6URNEO  /  LENGTH 


TIME  HEAD  EXIT 

P*  PI*  ISP 

PLfcN 

IMP  SUM  W  OEPL 

0.0 

1.00 

16.00 

20.00 

24.00 

28.00 

32.00 

36.00 

38.00 

40.00 

LOOP 

4.60  61.1  40.9 

40.9  71.9  1  50. 

14.0  1216.  8.057 

9  7948.  45.974  0.0 

2.084 

2.219 

2.325 

2.376 

2.421 

2.438 

2,438 

2.438 

2***3I 

2.292 

12 

4.70  61.9  31.0 

31.0  54.4  140. 

14.0  069.  6.196 

3  8052.  46.686  0.0 

2.111 

2.266 

2.351 

2.401 

2.438 

2.438 

2.438 

2.436 

2.438 

2.312 

12 

4.80  *.9,2  24.7 

24.7  43.3  129. 

14.0  646.  4.994 

3  8128.  47.246  0.0 

2.135 

2.270 

2.373 

2.422 

2.438 

2.438 

2.436 

2.436 

2.438 

2.328 

13 

4.90  40.2  20.2 

14.0  484.  4.122 

2.157 

2.291 

2.393 

2.438 

2.438 

2.431 

2.438 

2.431 

2.436 

2.339 

15 

35.5  117.5  8185.  47.702  0.0 


16 

- 1 .01166-02 

3.24456 

01 

1.0245E  00 

1.0003E 

00 

17 

-9.2380E-Q3 

5.24436 

01 

1  .02*56  00 

1.00006 

00 

\6 

-8 .43936-03 

3.24406 

r»l 

9.74526-01 

1.00006 

00 

19 

-7.75356-03 

3.24366 

01 

9.74526-01 

1.00006 

00 

20 

-7.1572E-03 

3.24376 

01 

9.74526-01 

1.00006 

00 

21 

-6.5607E-01 

3.24356 

01 

9 . 7452E-01 

l.OOOOE 

00 

22 

-6.0819E-03 

3.2434E 

01 

9.7452E-01 

l.OOOOE 

00 

23 

-5.76136-01 

3.2433E 

01 

9 ■ 7452E-01 

l.OOOOE 

00 

24 

-5.44266-03 

3.2432c 

01 

9  «  7  4  52  E— 0 1 

1.00006 

00 

25 

*5.1  270E-Q3 

3.24311 

01 

9.74526-01 

l.OOOOe 

00 

L«  WARNING 

5 .00  32.% 

16.4 

14 

0  343. 

3.353 

2. IT*  2.310  2.409  2. 43*  2*438  2.438  2.438  2.438  2.438  2.347 


16.4  21.7  102.3  8224.  48.075  0.0 

TMK0A7  UNCHCMfcO 


AVERAGE  THKUST  AND  ISP  ARE  1418.  LBM,  170.4  SEC  WITH  ACTION  TINE  OP  5.099  SECONDS 

MAXIMUM  THRUST  IS  2280.  LB.  ANO  MAXIMUM  PRESSURE  IS  941.42  PSU 


5.10  23.3  13.9  14.0  173.  2.348 

13.9  20.6  73.3  8252.  48.340  0.0 


2.192  2.325  2.419  2.438  2.438  2.438  2.438  2.438  2.438  2.354 

STOP  ISP  -  164.16,  BASEO  Oft  LOAOEO  HEIGHT  OP  50.3  L8S 


t 
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The  axial  distribution  along  the  port  at  t  =  0.10  shows  velocity  V 
increasing  to  sonic  (39,529  in. /sec)  at  node  17.  The  throat  is  beyond  the 
minimum  area  due  to  mass  addition,  and  area  A  at  node  17  is  not  the 
minimum  due  to  burning. 

The  run  proceeds  normally  with  pressure  decreasing  and  the 
throat  moving  aft  until  3.30  seconds.  At  this  time,  a  large  portion  of  the 
web  burns  out  and  the  L*  warning  is  bypassed.  Because  burnout  was 
sudden,  the  motor  pressure  decays  rapidly,  the  motor  unchokes,  and 
action  time  is  sensed  prior  to  the  3. 50- second  printout.  The  run  is  termi¬ 
nated  after  3.60  seconds  as  the  pressure  falls  within  2.0  psi  of  ambient. 

The  integral  of  mass  flow  and  the  loaded  weight  do  not  agree  well  for  this 
type  of  sharp  decay,  and  the  user  may  need  to  reduce  time  steps  during 
tailoff  for  more  precise  integration. 

The  second  case  in  the  sample  batch  has  a  cylindrical  port  and 
two  propellants;  the  aft  propellant  is  assumed  to  have  half  the  base  burn 
rate  of  the  forward  (and  prior  case)  propellant.  Also,  the  erosive  burning 
correlation  is  assumed  and  is  input  the  same  as  for  the  first  propellant. 

The  axial  length  of  the  beginning  of  the  second  propellant,  Z(6)  =  36,  "keys" 
the  option  and  is  noted  by  two  printouts  of  the  line  "EROSION  MODEL.  .  .  .  ". 

Since  other  changes  were  not  input,  the  second  case  is  the  same 
as  the  first  except  that  the  cylinder  of  propellant  No.  2  replaces  the  aft 
4  inches  of  propellant/cone  in  the  first  case.  The  slower-burning  aft 
propellant  causes  the  motor  to  tail  off  with  0.586  inch  of  web  remaining  and 
the  run  terminates  at  2.7  seconds  due  to  low  pressure. 

The  third  case  is  a  motor  with  fuse  ignition.  The  input  includes 
revised  geometry  to  set  the  port  ID  to  fuse  diameter  (table  3),  and  a  nega¬ 
tive  value  of  X2  in  Z(6)  is  input  to  cancel  the  two-propellant  option  from 
the  previous  case.  The  fuse  velocity,  VFZ  in  Z(15),  "keys"  the  fuse 
option  and  the  beginning  location  of  the  ignition  front  is  given  as  node  19 
in  NUM(9),  to  designate  starting  at  the  aft  end  of  the  grain.  The  time 
steps  for  this  ignition  option  are  internally  calculated  so  that  the  flame 
front  moves  one  node  for  each  step.  Because  the  web  burned  is  printed 
only  at  alternate  nodes,  the  flame  front  moves  one  print  increment  for 
every  other  time  step.  After  the  fuse  burns  to  the  head,  the  time  steps 
are  obtained  from  the  tabular  input  as  in  normal  runs. 

The  rise  time  is  input  as  0.8  second  in  Z(  2)  to  avoid  the  L* 
instability  test  during  the  fuse  ignition  period.  (A  warning  would  occur 
at  0.02  second  and  the  test  would  be  canceled,  so  the  meaningful  warning 
at  5.0  seconds  would  be  missed.)  The  trial  pressure  PSTART  is  input 
as  14.0  psia  in  Z(22),  and  with  this  option  becomes  the  initial  pressure 
in  the  port  for  a  degenerate  calculation  at  ignition  time. 


AFRPL-TR  -73-20 


The  printout  of  parameters  along  the  port  preceding  the  0.12- 
second  normal  output  indicates  that  the  flame  front  has  moved  six  incre¬ 
ments  along  the  motor,  consistent  with  the  web  burned  profile  for  that 
time. 


The  calculations  do  not  converge  at  the  0.20- second  point  but  are 
oscillating  between  808  and  813  psia.  The  computational  stability  is 
particularly  unpredictable  because  a  wide  range  of  geometry,  erosion 
model,  burn  rate  slope,  grain  deflection,  or  aluminum  fraction  burned 
are  considered.  In  this  example,  the  5-psia  uncertainty  may  be  assumed 
small  in  the  effect  on  overall  performance  parameters.  The  user,  how¬ 
ever,  must  observe  that  results  are  within  requirements.  At  5.0  seconds, 
the  calculation  is  substantially  overdamped,  so  decreasing  the  loop  gain 
FACTOR  in  Z(25)  will  reduce  the  oscillation  at  0.20  second,  but  will 
increase  the  number  of  iterations  later.  A  desirable  solution  is  to 
decrease  the  time  step  size  during  tailoff,  then  decrease  FACTOR. 

The  run  is  terminated  due  to  the  input  of  TSTOP  =  5.0  seconds 

in  Z(2). 


(The  reverse  is  blank) 
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APPENDIX  A 

SYMBOL  LIST  AND  SOURCE  PROGRAM  LISTING 


A  complete  list  of  symbols  employed  in  the  FORTRAN  program  is  given  in 
Table  A-l.  The  list  is  approximately  in  alphabetical  order  but  parameters 
are  grouped  to  illustrate  similarity  ol  related  terms.  For  example,  all  of 
the  symbols  representing  burn  rate  are  grouped  under  R.  Any  parameters 
that  are  input  are  identified  in  the  second  column  of  the  table. 

The  source  program  listings  of  the  main  program  and  subroutines  DATA  IN, 
TABIN,  and  TAB  comprise  Table  A-2.  The  user  may  find  it  convenient  to 
recognize  the  correspondence  between  the  EQUIVALENCE  statements  in  the 
main  program  and  Tables  IV-3  and  IV-4  in  Section  IV  of  the  User's  Manual. 


A-l 
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TABLE  A-  1 
SYMBOL  LIST 


Symbol 

Input 

Pi  rimrtrr 

Unit  • 

Dr  nr  r;ptlon 

A 

A 

in.2 

Arran,  port,  to  trap  slivers 

AL 

in  2 

■  avrd  A  last 

AP 

in.1 

instantaneous  port 

ABAR 

in. 

average  AP  (I)  and  AP  (J) 

AHEAD 

in  2 

hrad  r nd  of  g rain 

ADFFL 

in  2 

drflrctrd  port 

APT 

in.2 

■  avrd  throat  area,  oard  in  pressure  guess 

ACT 

- 

Code  word,  action  timr  printout  logic 

.. 

ALPHA 

Factor,  Lenoir  -  Robilla  rd  (  L- R )  rq 

X 

ALPHAS 

- 

Relative  thermal  diffusivity 

Al 

X 

ALUM 

- 

Aluminum  fraction  of  propellant 

ANG 

- 

Divergence  angle,  thrv&t  calculation 

ATERM 

Factor  in  L-R  eq,  oCi',#/Xfcr2 

BTERM 

- 

F  tetor  in  L-R  eq,  -  Pp/G 

l' 

X 

BETA 

- 

Constant  in  L-R  rq 

BE  TAR 

1 b/ in } 

'  P 

X 

Bll 

- 

First  coefficient,  L-R  eq 

d 

C 

in. /ire 

Sonic  velocity 

X 

CHOKE 

- 

Code  word  to  indicate  choking 

COE  EE 

Fdf!or,  shear  stress  calculation 

X 

CON  V 

Convergence  ratio 

CON  Vi 

- 

l  •  CONV,  axial  increment* 

CON  V  1 

- 

10  •  CONV,  choke  constraint 

COUNT 

Counter  for  table  lookup,  time  of  detailed  output 

COUNTP 

time  of  sensitivity  tabulations 

( 

P.  K 

X 

CPG 

kcal/g-mde 

Specific  heats,  gas 

c  c 

p,  P  or  p 

X 

C  PP 

kcal/g-mole 

products  (when  f  =  1.) 

C  1 

_ 

1  -  •  (1  -  f) 

o 

Cl 

- 

l  -  •  -  f1) 

X 

I) 

in. 

Diameters,  port  innt  r  input 

X 

DM  AX 

in. 

grain  outT 

DIAM 

in. 

instantaneous  port 

X 

DPLF.N 

in. 

plenum 

X 

DTH 

in. 

♦hroat  of  nozzle 

X 

DALUM 

mic  rons 

aluminum  particle 

DR 

in.  / sec 

Increments,  erosive  burning 

DP 

psi 

pressure 

DU 

in.  / sec 

veloc  ity 

DW  E  n 

in. 

web  burned 

DELL 

in. 

aft  end  burning  length 

DE  LP 

psi 

pressure  r  r  ror 

X 

nr  lt 

h  l*( 

time  int  rements 

1)1  LT  l 

R  Cl 

fct/2 

n / 

in. 

length 

DX 

- 

L/N 

DXDT 

in.  »et 

fc  x  fc  t 

A-2 
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TABLE  A-  1  (Continued) 


Sym  ►xj  1 

Input 

P*  rimftr  r 

Unit* 

Dr  sc  ription 

DAP 

in.1 

Incremental  area  change,  used  alao  in  locating  throat 

DD 

in. 

Diameter  of  burn  beyond  outer  diameter  (for  aliver  log! 

DDJ 

in. 

[,a*t  atation 

DEFL 

in. 

Deflection,  instantaneous 

X 

DEFLO 

in. 

reference  input 

OF  FLN 

in. 

transient  ratio 

DFNOM 

• 

Term  in  throat  location  logic 

DIVR 

- 

Divergence  loss  term  in  thrust  calculation 

X 

DOME 

- 

Key  to  hemispherical  head  option 

F 

_ 

Error  ratios,  saved  value  of  inner  loop 

ERRM 

- 

Mach  number,  outer  loop 

ERRML 

- 

last  value  ERRM 

ERRN 

- 

over-choke  ratio 

FRRU 

- 

velocity 

• 

EPS 

- 

Solid  fraction,  with  slip 

EPSI 

- 

1  -  ♦ 

’  o 

EO 

- 

Solid  fraction,  without  slip  (f  =  1 ) 

EOI 

- 

l  -  •  o 

X 

EPSP 

- 

Plenum-to-throat  expansion  term,  nozzlrd  motors  only 

EXPWNT 

- 

-  w  t 
e  n 

r 

F 

FI 

- 

Slip  ratio,  u_/u 

1  -  f  V  * 

FISQ 

- 

1  -  f1 

FAC 

. 

Convergence  factor,  outer  loop  gain 

X 

FACTOR 

- 

input  value 

FAULT 

- 

Indicator  of  choke  constraint  error 

F 

FG 

lb 

Thrust,  instantaneous 

FGL 

lb 

save  prior  value 

FAV 

lb 

average  over  action  time 

FM  AX 

lb 

maximum 

FAL 

Factor  in  aluminum-burned  calculation 

v 

FR 

- 

Fraction  of  aluminum  burned 

FILL 

lb 

Loaded  propellant  weight 

FIX 

- 

Term  used  to  gue  is  next  head  pressure 

Fx 

- 

- 

Wall  frictional  force  per  unit  length 

r. 

G 

lb/ sec  -  in.2 

Specific  mass  flow 

>. 

N  n 

. 

Ratio  of  specific  heats  for  sonic  velocity  calculation 

GAMP 

- 

Ratio  of  specific  heat  terms,  products 
=  Tc  -  Ts/u2  (energy  equation) 

r 

GAMT 

- 

GFI 

- 

V  -  1 

GF2 

- 

V  ♦  1 

GFI 

- 

v  -  1/ 

GFS 

- 

i/i  >  ♦  i) 

GO 

in.  sec 

Mass/weigi.t  conversion 

°o 

GOi 

lb/  in.2  -  sec 

Erosive  coefficient 

GOiX 

r  m 

°0 
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TABLE  A-  1  (Continued^ 


Symbol 

Input 

Parametc  r 

Input 

ttesi  ription 

H 

Matrix  tn  save  pararnrtric  results 

X 

HIGH 

in. 

Height  of  two-dimensional  t  ,r\ 

I 

Index,  axial  nod**  numbr r 

X 

1GN 

Ignition  time  option  key 

IPR 

Axial  print  index 

IS 

Index,  time  step  in  pa  ramet  r it  run 

I 

HP 

ISP 

R  Cl 

Specific  impulse 

J 

Index,  I  -  1 

k 

KI) 

lb/  in.' 

Spring  constant  in  transient  deflection  ralculatioi 

L 

Index,  instantaneous  head  end 

LOOP 

Ix)op  counter,  inner  loop 

LOO  PC 

outer  loop 

LOOP  L 

limit,  inner  loop 

LOO  PM 

limit,  outer  loop 

1. 

LSTAR 

in. 

Characte  ristic  length 

X 

MO  LG 

Moles  of  gas/  100  gr 

X 

MOLP 

Moles  of  products/ 100  gr 

Mu 

X 

MW  TG 

Molecular  weight  of  gas 

MOV  V 

Code  word  to  trap  throat  movement 

rh  ( 

. 

lb/in.  -sec 

Mass  flow  added  per  unit  length 

M 

MU 

lb  in. -  sec 

Viscosity  of  gas 

MVZ 

„0.l 

X 

N 

Number  of  axial  nodes 

X 

NXRAY 

Index,  location  of  r  output  «n  parametric  runs 

X 

NSTAR 

Index,  throat  location 

NP 

Index,  propellant  number 

NPAR 

Index,  parameter  number  in  parametric  runs 

P 

X 

P 

psia 

Pressure,  local 

PL 

psia 

last  value  of  p  .  . 

PAMB 

psia 

,  •  ,  n  e  a  cl 

ambient 

PLO 

psia 

last  value  at  head 

X 

PPLEN 

psia 

plenum,  only  if  noz/.lrd 

X 

PR  FF 

psia 

reference  for  deflection  calculation 

X 

PST  Ap  T 

psia 

t  rial  value 

PM  AX 

psia 

maximum  at  head 

Pr>00 

- 

p/SOO 

PT  STAR 

psia 

throat  stagnation 

PA 

- 

term  in  star  fuse  geometry 

P 

P  DOT 

psi,  set 

Rate  of  head  pressure  change 

PNUM 

Floated  NP,  propellant  number 

r 

PI 

. 

1. 14  159 

P  U 

- 

"/<? 

PI4 

- 

!»/  4 

PIfl 

- 

•v  H 

l'r 

PR 

- 

Prandtl  number 

A-4 
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TABLE  A- 1  (Continued) 


Input 

Pa  rametr  r 

Unit  • 

Drat,  nption 

r 

R 

in .  /ifc 

Burn  rate,  inatant jneoui 

RL 

in.  / »rc 

saved  R 

RREAL 

In,  /ire 

t  rapped  for  alive ra 

R  I 

in.  iff 

laat  v  »lue  saved  to  amooth  v  ill  shear  atreaa 

rri 

in.  /  nee 

present  value  to  smooth  wall  shear  stress 

r 

o 

RO 

in.  ire 

base  bur n  rate 

RATIO 

MOLCl  MOLP 

RHI 

in. 

Hydraulic  radius,  ith 

RHJ 

in. 

jth 

X 

KHO 

lb/  in.1 

Density,  aolid  propellant 

, 

RHOG 

lb/in.1 

R  *  • 

RIIOL 

lb  in.* 

gas,  saved  value 

KHOBAR 

1  b/  i  n  } 

average  of  ith  and  jth 

X 

RHOAMR 

lb/  in,1 

initial  condition  in  motor 

R 

RG 

in./°i 

Gaa  constant,  gas 

R 

RP 

in./°F 

products 

R 

u 

RU 

in. /UR 

universal 

RT 

in. 

Product  Rp  •  Tc 

S 

in. 

Perimeter  of  port 

X 

SALUM 

Coefficient,  aluminum  slip  correlation 

SAL 

- 

Term,  aluminum  slip  correlation 

fAl 

SUP 

- 

Aluminum  slip  ratio 

X 

SALOX 

. 

Coefficient,  AljOj  slip  correlation 

SOX 

- 

Term,  A120»  slip  correlation 

fAljO, 

Sl.IPOX 

Al|Oj  alip  ratio 

STAR 

- 

Number  of  points  in  star  point 

STEP 

- 

Counter  in  parametric  runs  -  =  omit,  0  =  base  run,  ♦  =  counter 

SUMr 

lb 

Running  sum,  thrust 

SUM  W 

lh 

mass  flow 

TAIL 

- 

Indicator  for  tailoff  logic 

TANN 

* 

Tan  (star ) 

▼ 

TAUW 

psi 

Wall  shear  stress, 

T 

TAUBAR 

pai 

average  in  Ax  increment 

tF 

X 

TC 

°R 

Chamber  temperature 

T 

rs 

°R 

Local  static  temperature 

t 

TIME 

see 

T  imi*, 

X 

TSTOP 

ser 

completion 

X 

TPART 

sec 

values  at  parametric  points 

X 

T  LIST 

sec 

valuta  for  detail  output 

TIME1 

aec 

adjusted  to  avoid  rounding 

TEMP 

Temporary  storage  of  factor  tc  be  used  with  several  statements 

u  or  U 

U 

in.  sec 

Gas  velocity,  instantaneous 

R  K 

UTRY 

in. / sec 

for  convergence  check 

UMAX 

in. / sec 

at  T  —  0  to  trap  supersonic  flow  logic 

up 

- 

Particle  velocity 

iJi'W 
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TABLE  A-l  (Continued) 


Symbol 

Input 

Parameter 

Units 

V 

in? 

VOL 

in? 

VOLL 

in? 

X 

VFY. 

i**,  /sec 

W 

lb/ sec 

WL 

lb/ sec 

WADD 

lb/ sec 

W  BAR 

lb/ sec 

w. 

X 

WIGN 

lb/ sec 

ign 

W,  7. 

lb/sec 

WNL 

lb/sec 

.’.'BP 

in. 

V.  EBM 

in. 

X 

WIDE 

in. 

wn 

X 

WN 

rad/sec 

WNT 

- 

WNT  1 

- 

WN2 

- 

WN2T 

- 

X 

X 

in. 

xo 

in. 

X 

X2 

in. 

XE 

in. 

XHEAD 

in. 

X 

XPLEN 

in. 

XL 

- 

XSO 

- 

XD 

in. 

XDO 

in. 

XDOT 

a 

X 

XKG 

in.1/ lb 

m 

XM 

- 

XR 

in. / sec 

XTERM1 

to  5 

Y 

in. 

ZTH 

i  . 

_ Pete  r  lpt  ton _ 

Volume,  running  sum  along  axis 
of  plenum 
caved  value 

Velocity  of  ignit ion  upet  ream  spreading 

Mail  flow,  instantaneous 
caved  value 

added  In  Ax  increment 
average  inAx  increment 
igniter  or  head  end  dome 
fuce  igniter 
caved  value,  exit 

W  eb  bu  rned 

Initial  or  maximum  web 

Initial  length  of  ctar  fine 

Natural  frequency,  deflection  calcul/  tion 
w,,  •  6  t 
1  -  wn  fc: 

2/  wn 

2/ wn  ♦  6 1 

Length, 

to  inctantaneouo  head 
to  oecond  propellant 
motor 

uninhibited  head  web  burn 

plenum 

X/L 

x/so 

Deflection  transient  term 
prior  value 

Rate  of  change  of  XD 

Coefficient  in  erosive  correlation 

Exponent  in  erosive  correlation 

Temporary  value  of  r 

Factors  of  transient  deflection  calculation  equation 
Instantaneous  wid'h  of  star  fin  port 
Width  of  two-dimensional  throat 
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TABLE  A-2 


SOURCE  LISTING 


u  o  o  o  o  •_> 


r 
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02/20/73  PAGE 


C 

C 

C 

C 

C 

C 


NQZZLELESS  RUGAE T 
DYNAMIC  MODEL 


MOTOR  BALLISTICS  /  EROSIVE  BURNING 


LOCKHEED  PROPULSION  COMPANY  /AIR  PURGE  ROCKET  PROPULSION  LAB 


MWTGt  MU,  MU2 , 


LSTAR 


REAL  ISP,  KO,  MOLG,  MOLP, 

DIMENSION  LOGIC  1 25 1 ,  NUMI25I 

COMMON  /COM2/ZI125)  /COM5/SVSTEM 1125  )  /COM1/COMI 102 ) 

I  /C0M4/HAS  T  ER ( 3000 ) 

EQUIVALENCE  ( LOGIC ,SYST EMI  101 »»  ,  (NUM, 2(101)1 

DIMENSION  X  (  30 1 ,  HEBPI30),  DMEBI30),  R(30),  RLI30),  PI  30  >  t 


1  U(  30  >  ,  A  (  30 ) ,  P.H0GI30),  GOO),  0(30), 

1  MEBMI3D,  TAUNI30),  RHOSARI 30 ) ,  M(30), 

1  RHOL  (301,  AL  (  30 ) ,  DEFLOOO),  DEFLI30), 


AP ( 30 ) ,  DMAX ( 30 ) , 
Cl  30  >  •  E I  30 )  , 
OIAMOO),  HI  9,  1 5  ) 


DIMENSION 
1  BETARI2  > , 


XKG( 2 ) , 
XM ( 2  I , 


BH(  2 )  • 
GO  2X12) 


GO 2 ( 2  )  ,  PR ( 2  )  ,  RH0I2),  ALPHA (21, 


(LIST ,NUM ( 12)1 


624-0002 
PROPELLANT  624-000  3 
624-0004 
624-0005 
624-0006 
624-0007 
624-0008 
624-0009 
624-0010 
624-0011 
624-0012 
624-0013 
624-0014 
624-0015 
624-0016 
624-0017 
624-0018 
624-0019 
624-0020 
624-0021 
624-0022 
624-0023 
624-0024 
624-002  j 
624-0026 


EQUIVALENCE  (N,NUN) , 

N  »  10,  19,  28,  ONLY 
NUM 1 4 )  *  -1  IS  A  REQUIRED  KEY 

( KOOE, NUM( 21),  IKEV,NUM( 5) I •  I  IGN,NUMI6)  )  ,  (NSTAR,NUM( 7  I ) 

1  (  NXRAV ,NUM 18)),  (L ,NUM( 9 ) ) ,  (LOOPM, NUMI  10 )  ) 

EQUIVALENCE  ( TSTOP, 2 )' »  ( TR1SE ,2(2  )  )  ,  (ALUM, 2(3)),  (BETA, Z(4) ) ,624-0027 

1  (MU, 2(5)),  ( XE ,2 ( 7 ) ) ,  (XPLEN,2(8) ) ,  IDPLEN,2I9) ) ,  624-0028 

1  I DTH, 2Th,2 ( 10 ) ) ,  (PPLEN  ,7(11>),  (MIOE.ZI  12) ),  ( HI GH, 2 ( i 3 ) ) , 

1  (STAR, 2(1411,  (VE2, 2(15)1,  (STEP, 2 ( 18  ) )  ,  (DOME, 2(19)) 

EQUIVALENCE  ( EPSP, 2  ( 1 7 ) )  ,  (  X2 , 2  6 ) I ,  ( CHOKE , Z ( 21 ))  , 

1  (FACTOR, 2(251),  (WFZ, 2(16)1 

EQUIVALENCE  (CONV, 2< 20 ) ) ,  ( RH0AM6 ,2(23)  )  ,  I PSTART, 2 ( 22 ) )  , 

1  (PAMB, 2(24)  I  ,  (KD, 2(261),  (WN,2I27)),  ( PREF,Z I  28 ) ) , 

1  ( AL  PHAN, 2(29) ) , 

l  ( G02 ,BH,Z (30)1,  I  XKG , PR , 2  I  32 ) ) ,  (XM,Z(34>),  IRH0,2(36»), 

1  ( OALUM, Z ( 38  )  ) ,  ( SAL  OX ,2(39) ) ,  ( S AL UM, 2 ( 40 )  ) 


EQUIVALENCE  (IS, 2(50)),  ( DELL , 2 ( 51  I  )  ,  (WL, 2(521), 


624-0029 
624-0030 
624-0031 
624-0032 
624-0033 
624-0034 
624-0035 
624-0036 
624-0037 
624-0038 
(  XDO, 2(53) ) ,624-0039 


( XOOT, 2 ( 54  |  )  , 

( PMAX , Z ( 58  )  ), 

(COUNT, 2(62  )  ),  (COUNTP, 2(631 ),  I T A IL , Z ( 64 ) ) , 

(ACT, 2(66)),  ( FGL ,Z ( 67) I  ,  ( XHEAU, 2(68  )  ) 


( FILL ,Z( 55  )  I  , 
(HNL, 2(59)1, 


(SUMF ,2 156) )  , 
(TIME, 2(60)  I, 


( FMAX, 2 (5711, 

( DEFLN, 2(61  )  >, 

(SUMW,ZI65)>, 


GO  - 
PI  ■ 
P 12  > 
P 1 4  « 
P  I  8  < 
RU  « 
00  7 


386.09 
3.141592 
PI  /  2. 

PI  /  4. 

PI  /  8. 
1545.  *  12, 
1*1, 3000 


7  MASTER! I ) 


DO  8 
2(1) 


1-1,125 

'  0. 


624-0040 

624-0041 

624-0042 

624-0043 

624-0044 

624-0045 

624-0046 

624-0047 

624-0048 

624-0049 

624-0050 

624-0051 

624-0052 

424-0053 

624-005* 
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00  9  1-1,102 


624-005* 


9  COM) ||  -  0. 

READ  NEW  HEADING  ANO  OAT  A  OR  CHANGES  TO  DATA 
UPDATE  INPUT  LOGIC  WORDS 

1  READ  (1,991  I 
IE  (I.LT.ll  CALL  EXIT 
WRITE  <3,991  I 

99  FORMAT  (  I  1 , 7 1 H 

1  I 

CALL  TAB  I N 
DO  2  1-1,100 

2  IF  (SVSTEMd  NE.O.I  2(11  -  SYSTEH(I) 

00  3  1-1,25 

3  IF  (LOGIC ( I  ) . NE .0 )  NUMI  I)  -  LOGIC ( 1 ) 

IF  ( NUMI 6  I .EO.O I  STEP  ■  -1. 

NXRAV  -  NUMI  8 1 

IF  (NXRAV, EO.O)  NXRAV  -  N  -  4 

4  00  5  1-50,68 

5  2(1)  -  0. 

CALL  TAB  (COUNTP,  0.,  TPART ,  11) 

TANN  -  TAN  (STAR  /  <2.*  Pill 
NS  TAR  -  NUM( 7) 

MOVE  -  NS  TAR 

IF  (NSTAP.EO.O)  NS  TAR  -  N 
L OOPM  •  -.AXO  INUMI10),  10) 

IGN  -  NUMI 6 ) 

IPR  -  N  /  9 

CONV  >  A MAXI  (CONV,  1 .6-5 ) 

C0NV2  -  2.  *  CONV 
CON VI  -  10.  *  CONV 
WF2  •  2(161 
CHOKE  -  2(211 
IF  (KD.EQ.O.)  KO  -  1. 

IF  (WN.GT.O.)  WN2  -  2.  /  WN 

C  SET  FOR  ROUND  OR  2-DIMENS IONAL  PLENUM  AND  THROAT  SECTIONS 

C 

AP(N*2)  «  P 1 4  *  DTH**2 
AP(N*l)  -  P 1 4  *  D  PLEN**2 
IF  (HIGH. EO.O.)  GO  TO  6 
APIN+21  -  HIGH  *  2TH 
API N+l I  «  HIGH  *  D  PLEN 

6  VOLL  -  AP(N*i)  *  X  PLEN 
VOL  «  VOLL 

PA  -  AMAX1  (STAR,  1.) 

C 

C  PREPARE  TERMS  FOR  TWO  PROPELLANTS 

C 

MU2  «  MU**. 2 

COEFF  «  .0288  *  HU2  /  GO 
X2  -  2(6) 


624-0056 

624-0057 

624-0056 

624-0059 

624-0060 

624-0061 

624-0062 

624-0063 

624-0064 

624-0065 

624-0066 

624-0067 

624-0068 

624-0069 

624-0070 

624-0071 

624-0072 

624-0073 

624-0074 

624-0075 

624-0076 

624-0077 

624-0078 

624-0079 

624-0080 

624-0081 

624-0082 

624-0083 

624-0084 

624-0085 

624-0086 

624-0087 

624-0088 

624-0089 

624-0090 

624-0091 

624-0092 

624-0093 

624-0094 

624-0095 

624-0096 

624-0097 

624-0098 

624-0099 

624-0100 

624-0101 

624-0102 

624-0103 

624-0104 

624-0105 

624-0106 

624-0107 
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NP  ■  2 


624-0108 


IF  (X2.GT.0.)  GO  TO  11 
NP  ■  1 

X2  •  X  E  ♦  1. 

11  00  13  1*1, NP 

IF  (XM< I ) . EO.O. I  XM(  11  -  1. 

G02XII)  *  &02II)**XM(1) 

BETARI I)  -  -  BETA  *  RHO< 1 1 

WRITE  13,121  NUN  I  I >1 1 ,  I,  RHO (II,  BHtl),  PR (I) ,  XMIII,  ZIH-5) 

12  FORMAT  (17H  EROSION  MODEL  NO  12,  15H  FOR  PROPELLANT  12,  10X, 

1  6H  RHO  *  F  7 .4 , 5X  *4HK 1  -  F7.4,5X,4HK2  •  F7 .4 , 5X , 5HEXP  «  F5.2, 

1  5X, 1 1HT0  LENGTH  -  F5.1.4H  IN.I 

13  ALPHAII)  -  .0288  *  MU2  /  PRI1I**. 6667  •  BH  (  I  1  /  RHOtll 
C 

C  FILL  LISTS  WITH  GEOMETRY  AND  TRIAL  VALUES 

C 

I  -  1 
NP  «  1 
XIII  -  0. 

TEMP  «  N  -  1 

CALL  TAB  10.,  2.,  OMAX 111,  31 
GO  TO  15 

14  XL  -  XI  1-1  1  /  XE 

CALL  TAB  I  XL ,  0.,  OX,  4) 

OX  *  OX  *  XE  /  TEMP 
XIII  -  XI 1-1  I  ♦  DX 
IF  (XIII.GT.X2I  NP  *  2 
CALL  TAB  (XI  I l,2.,0MAX( 11,31 

FILL  ■  FILL  ♦  PIB  *  OX  *  RHOINPI  *  (DMAX(I>**2  ♦  DMAX I l-l I **2 > 

15  CALL  TAB  (XI  1 1,1. ,01 11,31 
WEBMIII  *  .5  *  I OMAX I  1  I  -  0(1)1 
CALL  TAB  (XIII,  0.,  DEFLOI  1 1 ,  7) 

UIII  *  0. 

Ell)  *  0. 

Rill  -  0. 

RUN  -  0. 

WEBPI 11*0. 

GII)  »  l.E-6 
Will  *  0. 

PIN  «  0. 

API!)  »  0. 

TAUWll  I  *  0. 

RHOGII)  >0. 

I  «  1  +  1 

IF  (I.LE.N)  GO  TO  14 
OIAMINSTAR I  *  OINSTAR) 

PL  ■  P  PLEN 
PLO  *  PPLEN 

PIN+ll  «  A MAXI  IP  PLEN,  PAMB I 
AHEAD  *  P I  4  *  DMA  X(  1  1**2 
CALL  TAB  (0.,  0.,  OELT ,  5) 

CALL  TAB  (COUNT,  0.,  TL 1ST,  6) 


624-0109 

624-0110 

624-0111 

624-0112 

624-0113 

624-0114 

624-0115 

624-0116 

624-0117 

624-0118 

624-0119 

624-0120 

624-0121 

624-0122 

624-0123 

624-0124 

624-0125 

624-0126 

624-0127 

624-0128 

624-0129 

624-0130 

624-0131 

624-0132 

624-0133 

624-0134 

624-0135 

624-0136 

624-0137 

624-0138 

624-0139 

624-0140 

624-0141 

624-0142 

624-0143 

624-0144 

624-0145 

624-0146 

624-0147 

624-0148 

624-0149 

624HH50 

624-0151 

624-0152 

624-0153 

624-0154 

624-0155 

624-0156 

624-0157 

624-0158 

624-0159 

624-0160 
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WRITE  0.10)  X(NSTAR),  X(N),  DID.  DINSTAR).  DIN).  DMAX(l), 

1  UMAX (NSTAR ) ,  DMA X< N I ,  N 

10  FORMAT  I 18H0  NOMINAL  GEOMETRY/ 18X»4HHEA04X»6HTHR0AT6X,4HEX  IT/ 

1  8 HO  LENGTH1 4X.2F10 • 2/1 2H  INNER  DI AM3F10.2/ 12H  OUTER  DIAM3F10. 
1  47H0  DETAILED  GRAIN  DESIGN  INPUT  IN  TABLE  3  TO  USE  13, 

1  33H  NODES  SPACEO  AS  INPUT  IN  TABLE  4/1H0) 

WRITE  13,16)  NUMI5),  NXRAY,  IGN,  HIGH,  WIDE,  VF2,  DOME,  PPLEN, 

1  DALUM,  WN 

16  FORMAT  (42H0  CONTROL  WORDS  IN  INPUT  (NON-ZERO  IS  KEY)/ 

1  10X,  25H  PARAMETRIC  VARIATION  KEY  14,  5X,  13H  XRAY  NODE  NO  13/ 

1  10X,  25H  IGNITION  INPUT  DT  OPTION  14,/ 

1  10X,  24H  2-DIMENSIONAL  GRAIN  KEY  F7.3,  7H  (HIGH)/ 

1  10X,  24H  N-PUINT  STAR  GRAIN  KEY  F7.3,  7H  (WIDE)/ 

1  10X,  24H  FUSE  IGNITION  CONTROL  F7.1,  7H  (RATE)/ 

1  10X,  24H  HEAD-END  DOME  OESIGN  F7.1/ 

l  10X,  24H  PLENUM-THROAT  CALC,  P  *  F7.1/ 

1  10X,  24H  FRACTION  BURNEO  KEY  D  *  F7.1,  SH  MICRONS/ 

1  10X,  24H  GRAIN  DEFLECTION  KEY  WN  F7.1  ) 

WRITE  (3,17)  IX(I),I*1,N,1PR) 

17  FORMAT  ( 1H1 1 5X, 8HPRESSURE  8X, 6HTHRUS T3X,4HFL0W36X , 19HWEB  BURNED 
1LENGTH/2X.4HT IME4X ,4HHEAD4X ,4HEX I T4X  »4HPL EN25X , 10F7.2»2X,4HLD0P/ 
17X»2HP*6X»  3HPT*4X, 3HI SP5X, 3HIMP4X, 5HSUM  W4X.4HDEFL/) 

CALCULATE  COEFFICIENTS  FOR  PARTICLE  LAG  CURVE  FITS 

FR  -  1. 

SLIP  «  1. 

SLIPOX  -  1. 

IF  (OALUN.LE.O.)  GO  TO  18 
X50  -  XE  /  50. 

FAL  ■  (1.  -  .0084  *  OALUM)  *  X50**.6633 

SAL  ■  SALUM  *  X50**.1759  /  (OALUM  /  100.)**. 0842 

SOX  *  SALOX  *  X 50**. 02 34 

FILL  INITIAL  CONDITIONS  FOR  FUSE  CONFIGURATION 

18  L  -  1 

IF  (VFZ.LE.O.)  GO  TO  19 
L  =  NUMI9) 

IF  (L.EQ.N)  TIME  *  -  CONV 
OELT  »  ( X  ( L )  -  X(L-l))  /  VF7. 

19  P(L>  -  P  START 
FIX  *  2. 

SET  CONSTANTS  TO  BEGIN  AT  THE  INSTANTANEOUS  HEAD  END 

21  CALL  TAB  (TIME,  0.,  WIGN,  2) 

WIGN  >  A MAXI  (WFZ,  WIGN) 

FACTOR  a  2(25) 

IF  (TIME.GT.TRISE)  FACTOR  =  FACTOR  /  FR 
IF  (WN.LE.O.)  GO  TO  22 
WNT  a  WN  *  OELT 
WNT 1  a  l.  -  WNT 


624-0161 
624-0162 
624-0163 
2/624-0164 
624-0165 
624-01 66 
624-01 37 
624-01 66 
624-0169 
624-0170 
624-0171 
624-0172 
624-0173 
624-0174 
624-0175 
624-0176 
624-0177 
624-0178 
624-0179 
/  624-0180 
624-0181 
624-0182 
624-0183 
624-0184 
624-0185 
624-0186 
624-0187 
624-0188 
624-0189 
624-0190 
624-0191 
624-0192 
624-0193 
624-0194 
624-0195 
624-0196 
624-0197 
624-0198 
624-0199 
624-0200 
624-0201 
624-0202 
624-0203 
624-0204 
624-0205 
624-0206 
624-0207 
624-0208 
624-0209 
624-0210 
624-0211 
624-0212 
624-0213 
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EXPWNT  *  EXP  (-WNT  ) 

WN2T  *  WN2  ♦  OELT 

XTERM1  *  KD  *  (2.  *  WN  *  XOO  ♦  XOOT ) 

XTERM2  •  KO  *  WNT 1  *  XOO  -  PLO  *  (1 .  ♦  WNT  ) 

XTERM4  «  OELT  -  NN2 

XTERH5  »  (1.  -  WEBP(MSTAR)  /  WEBM(NSTAR) )  /  PREF 
22  BELT2  -  .5  *  OELT 
I  -  L 
J  *  I 
V  -  0. 

NP  ■  1 
OX  =»  0. 

PNOM  «  1. 

ERRN  -  1. 

FAULT  -  1. 

TAU  BAR  *  0. 

ERRML  »  ERRM 

PREPARE  PARTICLE  LAG  TERMS 

IF  (DALUM.LE.O.)  GO  TO  24 
P500  -  PIL)  /  500. 

FR  =  AM  INI  (FAL  *  P500*» . 1 151 «  l.) 

SLIP  a  AMI Ni  (SAL  *  P500**.2271  ,  1.) 

SLIPOX  •  AMI NI  (SOX  /  P500**.0381 ,  I.) 

UPDATE  THERMOCHEMICAL  COEFFICIENTS 

INPUT  TABLES  ARE  ASSUMEO  IN  TERMS  AND  UNITS  OF  ISP  PROGRAM 


CALL 

TAB 

(FR, 

0.  , 

TC,  8) 

TC  « 

US 

*  TC 

CALL 

TAB 

(FR, 

0., 

CPP,  9) 

CALL 

TAB 

<fr. 

0., 

CPG,  10) 

CALL 

TAB 

(FR, 

0., 

MOLG,  13) 

CALL 

TAB 

(FR, 

0., 

MOLP,  14) 

RATIO  >  l 

MOLG 

/  MOL  P 

CALL 

TAB 

IFR  , 

0., 

MWTG,  15) 

RG  ■  RU  /  MWTG 

EOI  ■  MWTG  *  MOLG  /  100. 

EO  -  I.  -  EOI 

F  *  SLIPOX  ♦  (I.  -  FR)  *  ALUM  /  EO  *  (SLIP  -  SLIPOX) 

FI  -  1.  -  F 

FI SO  ■  1.  -  F  •  F 

GAMP  *  CPP  /  (CPP  -  1.907  *  RATIO) 

EPS  *  EO  /  (F  *  EOI  +  EO) 

Cl  *  l.  -  EO  *  FI 
C2  *  l.  -  EPS  *  F1S0 
EPS!  «*  1.  -  EPS 
RP  *  EPS1  *  RG 

TEMP  *  CPG  ♦  (CPP  /  RATIO  -  CPG)  /  F 
C  CONS  »  TEMP  /  (TEMP  -  1.987)  *  RP  *  GO 
GFl  »  GAMP  -  1. 

GF2  «  GAMP  ♦  1. 


624-0214 
624-0215 
624-0216 
624-0217 
624-0218 
624-0219 
624-0220 
624-0221 
624-0222 
624-0223 
624-0224 
624-0225 
624-0226 
624-0227 
624-0228 
624-0229 
624-0230 
624-0231 
624-0232 
624-0233 
624-0234 
624-0235 
624-0236 
624-0237 
624-0238 
624-0239 
624-0240 
624-0241 
624-0242 
624-0243 
624-0244 
624-0245 
624-0246 
624-0247 
624-0248 
624-0249 
624-0250 
624-0251 
624-0252 
624-0253 
624-0254 
624-0255 
624-0256 
624-0257 
624-0258 
624-0259 
62*- -0260 
624-0261 
6t 4-0262 
624-0263 
624-0264 
624-0265 
624-0266 


non 
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GF3  *  GAMP  /  GF1 
GF5  -  2.  /  GF2 

CIN+2)  =  SORT  IC  CONS  *  TC  *  GF5) 

R  T  =  TC  *  R  P 

GAMT  *  .5  /  <GF3  *  GO  *  RG>  *  (1.  -  EO  *  F1S0)  /  E01 
UMAX  *  SORT  I  TC  /  GAMT) 

IF  INN.LE.O.)  GU  TO  30 
POOT  *  IPIL)  -  PLO)  /  OELT 

XTERM3  «  XT6RM2  ♦  XTERM1  *  DELT  ♦  POOT  *  WN2T 
XD  «  EXPWNT  *  XTERM3  ♦  PLO  ♦  POOT  *  XTERM4 
DEFLN  -  XO  *  XTERM5 
GO  TO  30 

trial  values  for  first  try  down  the  port 

25  OOJ  «  00 

RJ  ■  RRI 
RHJ  =  RHI 

ox  *  xm  -  xij) 

DXOT  =  OX  /  DELT 
IF  (U(I).GT.O.I  GO  TO  30 
G ( I )  ■  G(J) 

pm  -  pij) 
um  =  uiji 
RID  «  RIJ) 

30  LOOP  -  -2 

LOOPL  «  L00PM 

IF  II.GE.NSTAR)  LOOPL  -  2  *  L00PM 
IF  IX(]).LE.X2I  GO  TO  32 
NP  *  2 
PNUM  *  2. 

C 

C  LOCAL  PURT  AREA  AND  BURN  PERIMETER  CALCULATIONS 

C  BRANCH  FOR  OIFFEREN;  GEOMETRY  ASSUMPTIONS 

C 

32  IF  I IGN.GT . 0)  RL (I)  -  R(I) 

PWEBI1)  »  OELT  *  ( R(  I  '  ♦  RLIlll 
01  AMI  *  0(1)  +  OMEBI  I 

00  ^  01  AM  I  -  UMAX! I ) 

01  AMI  I  )  -  AMIN1  (  01  AMI  *  DMAXIIN 
XO  »  XII)  -  X (L  ) 

IF  IHIGH.GT.O.)  GO  TO  33 
IF  (MIOE)  34,34,35 

33  API  I )  ■  01  AMI  I )  *  HIGH 

S  =  2.  *  I  0 1 AMI  I )  +  HIGH  I 
GO  TO  36 

34  API!)  «  P 14  *  DIAMU)**2 
S  «  PI  ♦  01  AMI  1 1 

GO  TO  36 

35  Y  »  .5  *  01  AMI  I )  /  TANN 
W2  *  WIDE  -  Y 

S  =  2.  *  W2  ♦  P12  *  0 1  AMI  I) 

APII)  *  OIAMII)  *  I P 1 8  *  0 1 AM ( I )  ♦  W2  ♦  .5  *  Y) 


624-0267 
624-0260 
624-0269 
624-0270 
624—02  7 1 
624-0272 
624-0273 
624-0274 
624-0275 
624-0276 
024-0277 
624-0278 
624-0279 
624-0280 
624-0281 
624-0282 
624-0283 
624-0284 
624-0285 
624-0286 
624-0287 
624-0288 
624-0289 
624-0290 
624-029 L 
624-0292 
624-0293 
624-0294 
624-0295 
624-0296 
624-0297 
624-0298 
624-0299 
624-0300 
624-0301 
624-0302 
624-0303 
624-0304 
624-0305 
624-0306 
624-0307 
624-0308 
624-0309 
624-0310 
624-0311 
624-0312 
624-0313 
624-0314 
624-0315 
624-0316 
624-0317 
624-0318 
624-0319 
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C 

C 

C 

C 


36  RHI  *  AP I  I  )  /  S 

AttAR  *  .5  *  (APUI  ♦  API  J  )  I 
All)  >  A  BAR 

IF  ITIME.EQ.O.)  GO  TO  80 

PROPELLANT  BURN  RATE  SECTION 

EROSIVE  BURN  RATE  CODE  FROM  NUM12)  OR  NUMI3) 

CALL  TAB  (Pill*  PNUM,RO  » 1  I 
IF  (NUM(NP-U)  -  21  41,42.43 

41  DR  -  IGI I >**XMINP)  -  G02XINP) )  *  SORT  IUI1)  /  CIN+2)) 

GO  TO  44 

42  DR  ■  Gl  1  )  -  G02INP) 

44  RID  =  RO  ♦  XKG(NP)  *  A  MAXI  (DR,  0.) 

GO  TO  50 

43  IF  ll.EQ.L)  GO  TO  49 

A  TERM  =  ALPHAINP)  *  Gil)**. 8  /  X0«*.2 
B  TERM  -  -  BETA  *  RHOINP)  /  GII) 

RII)  -  A MAXI  (RII),  RO) 

46  XR  «  RO  ♦  A  TERM  •  EXP  IB  TERM  /  RII)) 

ERR  *  ABS  II.  -  XR  /  Rill) 

RII  )  «  XR 

IF  (ERR  -  CONV!  50,50,46 

PERFORM  SUMMATIONS  AT  EACH  AXIAL  NODE 

EXECUTE  TRAPS  TO  SMOOTH  THE  SL 1 VER ING-OFF  TRANSIENTS 

49  RII)  >  RO 

50  R  REAL  *  Rill 

IF  I  DO.  GT.0.0 )  R  REAL  «  .25  *  A  MAXI  IDIAMID  -  DU),  0.)  /  DELT2 
TS  «  TC  -  GAMT  *  Ul I ) **2 
IF  ( 1 .GT .NS TAR)  GO  TO  51 
RHOG(l)  -  P I  I  1  /  (RP  *  TS) 

GO  TO  52 

51  RHOGII)  -  GII)  /  Ull)  /  EPS1 

52  Cl  I )  «  SORT  <C  CONS  *  TS) 

U  TRY  «  0(1) 

IF  (I.EO.L)  GO  TO  57 
DAP  ■  API  I  )  -  API J ) 

IF  IDO.LT.O.)  GO  TO  54 
IF  (ODJ.GT.O.)  GO  TO  56 
TEMP  *  DO  /  (OD  -  DOJI 
GO  TO  55 

54  IF  (DDJ.LT.O.)  GO  TO  56 
TEMP  »  DOJ  /  (DD  -  DOJ) 

45  All)  a  Ad)  ♦  .5  *  DAP  *  TEMP 

56  W  ADD  ■  RHOINP)  *  (All)  -  ALII))  *  DXDT 
WII)  »  HI J)  ♦  H  ADD 

RMOBARII)  ■  .5  *  (RHOGII)  +  RHOG(J)) 

IF  ( IGN.LE.O )  W( I )  =  HID  -  I RH08AR ( I l-RHOL I  I  ) )  *  ABAR  *  DXOT 
!  F  (WII) .GT.O.)  GU  TO  58 
ML)  «  .75  *  PIL) 

SO  TO  22 


624-0320 

624-0321 

624-0322 

624-0323 

624-0324 

624-0325 

624-0326 

624-0327 

624-0328 

624-0329 

624-0330 

624-0331 

624-0332 

624-0333 

624-0334 

624-0335 

624-0336 

624-0337 

624-0338 

624-0339 

624-0340 

624-0341 

624-0342 

624-0343 

624-0344 

624-0345 

624-0346 

624-0347 

624-0348 

624-0349 

624-0350 

624-0351 

624-0352 

624-0353 

624-0354 

624-0355 

624-0356 

624-0357 

624-0358 

624-0359 

624-0360 

624-0361 

624-0362 

624-0363 

624-0364 

624-0365 

624-0366 

624-0367 

624-0368 

624-0369 

624-0370 

624-0371 

624-0372 
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C 

C 

C 


57 


58 


59 


61 


C 

C 

C 


62 


63 


TRAP  TO  AVOID  CROSSING  MACH  ■  1 
COMPLETE  CALCULATIONS  AT  THE  NODE 

IF  (OOME.LT.O.I  M I  CM  ■  RHOINPI  •  A  REAL  *  P12  *  DIAM|I)**2  ♦ 
RRI  -  R  REAL 
MIL)  ■  NIGN 

CALL  TAB  (TIME,  0.,  HEADW,  16) 

IF  (HEADW. EQ.O.)  GO  TO  58 
NIL )  ■  MIL  I  ♦  RHOINPI  «  RREAL 
XIL )  -  XHEAO  ♦  (RIL)  ♦  RLILI) 

OEFLII)  »  OEFLN  *  0EFL01I) 


•  (AHEAO  -  APIL )  )  * 
«  HEADW  •  DELT2 


HE  AON 


OEFLII)  /  D1AMI I ) ) 


ADEFL  ■  API  1 1  *  (1,  ♦  A,  * 

GII)  -  Win  /  ADEFL  «  E01 
TAUWIII  ■  0. 

IF  (I.GT.NSTAR)  GO  TO  59 
UC  >  GII)  /  RHOGII)  /  EPS1 
UII)  -  AMIN1  ICC),  UC) 

IF  (t.EO.L)  GO  To  62 
U  BAR  -  .5  *  (UII)  ♦  UIJI) 

ERRN  -  UII )  /  UC 

RRI  ■  (2.  *  RJ  ♦  R  REAL)  /  3. 

IF  (GII).LT.I.)  GO  TO  61 
TEMP  .  UII)  *  EXPIBETARINP)  *  RRI  /  G(  I )  I 
TAUWIII  •  COEFF  /  XO  *  (GII)  *  X0>**.8  *  TEMP 
TAU  BAR  -  (TAUWIII  *  TAUWIJI)  /  (RHI  ♦  RHJ I  *  DX 
DU  -  UII)  -  UIJI 

DP  «  UBAR  /  GO  *  ( RHOBARI  I )  *  DU  *  C2  ♦  WADD  *  Cl  /  ABAR)  ♦ 
Pill  »  A  MAXI  (PIJI  -  OP,  .7  *  PID) 

IF  II.LE.NSTARI  GO  TO  62 

DELP  ■  PID  -  RHOGII)  *  RP  *  TS 

UC  -  UII)  ♦  GO  /  GINSTAR)  •  DELP 

UC  -  AMIN1  (UII)  ♦  .25  *  (UMAX  -  UII)),  UC ) 

UII)  ■  A MAXI  (CINSTARI,  UC I 
ERRN  «  UC  /  UII) 

CHECK  CONVERGENCE  AT  EACH  AXIAL  STATION 


LOOP  »  LOOP  ♦  1 
IF  (LOOP. LT. 01  GO  TO  32 
Ell)  *  0. 

IF  (I.EQ.L)  GO  TO  68 
IF  (UII).LE.O.)  GO  TO  66 
ERR  U  -  UII )/U  TRY  -  1. 

IF  IABSI ERR  U) .LT .C0NV2 )  GO 
IF  (LOOP.LT.LOOPL )  GO  TO  32 
Ell)  ■  ERRU 


TO  63 


VERIFY  THAT  THE  THROAT  HAS  NOT  MOVED  INTO  A  DIFFERENT  STATION 


FAULT  «  FAULT  *  ERRN 
IF  (I.GT.NSTAR)  GO  TO 
V  »  V  ♦  ABAR  *  DX 


64 


624-0373 
624-0374 
624-0375 
624-0376 
WIGN624-0377 
624-0378 
624-0379 
624-0380 
624-0381 
624-0382 
624-0363 
624-0384 
624-0385 
624-0386 
624-0387 
624-0388 
624-0389 
624-0390 
624-0391 
624-0392 
624-0393 
624-0394 
624-0395 
624-0396 
624-0397 
624-039 B 
624-0399 
TAUBAR624-0400 
624-0401 
624-0402 
624-0403 
624-0404 
624-0405 
624-0406 
624-0407 
624-0408 
624-0409 
624-0410 
624-0411 
624-0412 
624-0413 
6 ?4-0414 
624-0415 
624-0416 
624-0417 
624-0418 
624-0419 
624-0420 
624-0421 
624-04 22 
624-0423 
624-0424 
624-0425 
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IF  (CHOKE. L T .0. )  GO  TO  68 
IF  (J.LE.HOVE)  GO  10  67 

DEMON  =  MAUD  /  M(I)  •  GF2  ♦  4.  •  IAU6AK  /  (Pill  ♦  P(J)) 

IF  (DAP  /  A  bAK  -  DENOM.LE.O.)  GO  TO  67 
NSTAR  ■  J 
I  »  J 

WRITE  (3,66)  NS  TAR ,  XI  I  I  ,  OAP 

>  FORMAT  ( 2  5HO  THROAT  LOCATION  AT  NODE  13,  10H  AT  LENGTH  F5.it 
1  13H  IN,  SLOPE  -  F6.2/I 
GO  TO  *9 

)  IF  (FAULi.LT. L,»  GO  TO  65 
IF  ( l.LT.N)  GO  TO  66 
GO  TO  83 
»  MOVE  *  NS  TAR 
NSTAH  -  N 
LOO  PC  -  0 
GO  TO  22 

I  IF  ( 1 .EQ.NSTAR  I  GO  TO  74 
I  I  ■  I  ♦  1 
)  J  «  I  -  1 

IF  (  I.Lfc.N)  GO  TO  25 

DETERMINE  PLENUM  COMO  IT  IONS  I  IF  ANY  I 

CALCULATE  CONVERGENCE  FACTORS  FOR  CHOKED  EXIT  TO  GRAIN 

IF  (L.EO.N)  GO  TO  81 
IF  (P  PLEN)  74,74,70 
I  CALL  TAB  ( P( N>1 ) ,  P  NUM,  RIN+l),  3) 

DELL  -  DELT  *  RIN+1) 

VOL  ■  VOLL  ♦  AP(N-H)  *  DELL 
EPSP  «  GF5**I1.  /  GF1I 
IF  (2(171. GT. 0.1  EPSP  -  Z(17) 

TERM  ■  DELT2  *  C(N*2>  *  A(N*2>  *  EPSP  /  VOL 

MIN-M  >  *  MIN)  *  PA  ♦  RHO(NP)  *  RIN+l >  *  (APIN+l)  -  PA  *  AP(NI) 

M  BAR  =»  .5  *  (M(N4l>  ♦  ML) 

IF  (TERM. GT. 0.9)  GO  TO  72 

PIN+1)  -  (PL  *  (1.  -  TERM)  ♦  DELT  *  RT  /  VOL  *  MBAR)  /  II.  ♦  TERM 
GO  TO  74 

PIN41I  -  MBAR  *  RT  /  (EPSP  *  APIN+2)  *  C(N*2)> 

IF  (FAULT  -  1.)  75,76,81 
ERR  M  *  1. /FAULT  -  1. 

FAC  ■  3.  *  FACTOR  *  AMIN1  (SORT  IERRM  /  .003),  .5  *  FAULT) 

GO  TO  78 

ERRM  -  P  ( N-f  1  )  /  P(NSTAR)  -  1. 

IF  (CHOKE.LT. 0.)  GO  TO  77 
IF  (ERRM.LE.O.)  GO  TO  73 
CHOKE  *  -1. 

NSTAR  «  N 
WRITE  (3,82) 

FORMAT  (17H0  THROAT  UNCHOKEO  /) 

GO  TO  77 


624-0426 
624-0427 
624-0428 
624-0429 
624-0430 
624-0431 
624-0432 
624-0433 
624-0434 
624-0435 
624-0436 
624-0437 
624-0438 
624-0439 
624-0440 
624-0441 
624-0442 
624-0443 
624-0444 
624-0445 
624-0446 
624-0447 
624-0448 
624-0449 
624-0450 
624-C451 
624-0452 
624-0453 
624-0454 
624-0455 
624-0456 
624-0457 
624-0458 
624-0459 
624-0460 
624-0461 
624-0462 
624-0463 
)  624-0464 
624-0465 
624-0466 
624-0467 
624-0468 
624-0469 
624-0470 
624-0471 
624-0472 
624-0473 
624-0474 
624-0475 
624-0476 
624-0477 
624-0478 
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73 


77 

76 


100 


79 


80 


C 

C 

C 


ERRM  «  U(I)  /  Cm  -  1. 

FAC  »  AMIN1  (-FACTOR  *  ERRM,  .81 
GO  TO  78 

FAC  ■  FACTOR  *  A8S(1.01  -  U(N)  /  C  <  N )  > 

DEL P  «  AMAX1  (-.5*  ERRM  «  FAC  I 
P(L>  -  P(L>  *  (1.  ♦  DEL  PI 
LOOPC  «  LOOPC  ♦  1 
IP  (LOOPC.LT. 01  CO  TO  22 
IF  (A8StERRMI.LT .C0NV1 I  GO  TO  79 
IF  (LOOPC. LT.LOOPMI  GO  TO  22 

WRITE  ( 3,100 1  LOOPC  ,  ERRM,  P(L),  FACTOR,  FAULT 
FORMAT  (14,  1 P6E12.4 I 

IF  (LOOPC .EQ.lOOPM+1 I  FACTOR  -  FACTOR  -  SIGN  1.05,  ERRM  *  ERRMLI 

IF  (L  I0PC.LT.2  *  L OOPM I  GO  TO  22 

IF  ( ABS  ( DEL  P )  .t.  I  .CONW  I  60  TO  79 

IF  <ABS(0ELPI.CT..05I  GO  TO  1 

IF  (I.GE.NI  GO  TO  81 

FAULT  -  1. 

GO  TO  48 

AMO  BAR) I  I  ■  RHO  A MB 

FILL  *  FILL  -  A8AR  *  OX  *  RHO(NP) 

GO  TO  79 

CHECK  FOR  CHOKING  AT  END  OF  GRAIN 


81 


C 

C 

C 


83 


84 


85 


86 


IF  (CHOKE. GT.O.)  GO  TO  83 
V  «  V  ♦  VOL 

IF  (FAULT.EU.l.l  GO  TO  83 
CHOKE  ■  1. 

GO  TO  65 

NORMAL  OUTPUT  AND  ADVANCE  TIME  INCREMENT S 

TIME1  *  TIME  ♦  l.E-4 
IF  (TIME1.LT. TUSTI  GO  TO  85 
COUNT  «  COUNT  +1. 

CALL  TAB  (COUNT,  0.,  TL 1ST ,  6) 

WRITE  (3,84)  (IT  TAUW(I),  G(  I  )  ,  U(I>,  P(!>,  APII),  W  <  I  )  ,  R(I), 
l  RHOG(I),  1*1 ,N I 

FORMAT  ( 1H01 5X, 3HTAU9X, 1HG1 1 X , 1 HV1 1 X, IHPl 1X,1HA1LX, 1HW 1 IX , 1HR 1 IX, 
L  3HRH0/  ( I10.1P8E12.4)  ) 

WRITE  (3,17)  (X(  I), l*l,N,IPR) 

00  86  I  ■  L  y N 

IF  (EID.NE.O.)  WRITE  (3,100)  I,  Ed) 

RLtll  >  Rill 
ALd)  *  Ad) 

0(1)  -  D'll  ♦  DMEB(I) 

RHOL(I)  *  RHOBAR ( I  ) 

WEBP(I)  *  AMIN1  (WEBP(l)  +  .5  *  DWEB(I),  WEBMd)) 

IF  (T1ME.LE.0.)  GO  TO  94 
IGN  *  0 


624-0479 

624-0480 

624-0481 

624-0482 

624-0483 

624-0484 

624-0485 

624-0486 

624-0487 

624-0488 

624-0489 

624-0490 

624-0491 

624-0492 

624-0493 

624-0494 

624-0495 

624-0496 

624-0497 

624-0498 

624-0499 

624-0500 

624-0501 

624-0502 

624-0503 

624-0504 

624-0505 

624-0506 

624-0507 

624-0508 

624-0509 

624-0510 

624-0511 

624-0512 

624-0513 

624-0514 

624-0515 

624-0516 

624-0517 

624-0518 

624-0519 

624-0520 

624-0521 

624-0522 

624-0523 

624-0524 

624-0525 

624-0526 

624-0527 

624-0528 

624-0529 

624-0530 

624-0531 


A  - 1  7 


o  o  o  n  non 
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thrust  calculat ion 

OUTPUT  SECTION 
DIVR  ■  1. 

ANG  *  .5  •  ( 0  I  AN  ( N  I  -  DIAM(NSTAR)) 

IF  (ANG. Lb. 0. )  GU  TO  120 
TbMP  «  XIN)  -  X(NSTAR) 

DIVR  *  .5  •  II.  ♦  TENP  /  SORT  (ANG**2 
FG  *  WIN)  *  U ( N I  /  GO  *  01 VR  *  Cl  ♦  A 
FNAX  *  A MAXI  < FG,  FMAX) 

PMA X  *  AMAX1  (P(L )  ,  PHAX) 

SUMMARY  TERMS  ANO  ACTION  TIME 


624-0532 

624-0533 

624-0534 

624-0535 

624-0536 

624-0537 

624-0538 

2  ♦  TEMP**2I)  624-0539 

ADEFL  *  (PIN)  -  PAMB)  624-0540 

624-0541 
624-C->4 2 
624-0543 
624-0544 
624-0545 
624-0546 
624-0547 
624-0548 

(UINSTAR)/C(NSTAR) ) **2  >  **GF3  624-0549 


ISP  *  FG  /  WIN)  624-0546 

SUMF  *  SUNF  ♦  DEL  T2  *  (  FG  ♦  FGL  )  624-0547 

SUMW  *  SUNK  ♦  DEL  T 2  •  (W(N)  ♦  WNL  >  624-0548 

PTSTAR  =  P(NSTAR)  •  (1.  ♦  .5  *  GF1  •  (UINST AR )/C ( NSTAR )) **2  )  **GF3  624-0549 

IF  (ACT.GT.O.)  GO  TO  125  624-0550 

IF  (FG.GT.0.1  *  FMAX)  GO  TO  125  624-0551 

FAV  *  SUMF  /  TIME  624-0552 

TEMP  *  SUMF  /  SUMW  624-0553 

WRITE  (3,123)  FAV,  TEMP,  TIME,  FMAX,  PMAX  624-0554 

123  FORMAT  (28H0  AVERAGE  THRUST  ANO  ISP  ARE  F6.0,5H  IBM,  F6.1,  4H  SEC  624-0555 
1  7X»  20H  WITH  ACTION  TIME  OF  F6.3,  8H  SECONDS/  5X,19H  MAXIMUM  THRU624-0556 


1ST  IS  F6.0,  2BH  LB,  ANO  MAXIMUM  PRESSURE  IS  F8.2.5H  PSIA/I  624-0557 

ACT  =  1.  624-0558 

624-0559 

PRINT  L»  WARNING  -  CANCEL  DURING  TAILOFF  624-0560 

624-0561 

125  IF  (TAIL.LT.O.)  GO  TO  130  624-0562 

IF  (TIME.LT.TRISEI  GO  TO  130  624-0563 

IF  (WEBP(L).LT.WEBM(L >)  GO  TO  128  624-0564 

WRITE  (3,127)  624-0565 

127  FORMAT  ( 38H0  L*  WARNING  SUPPRESSED  DURING  TAILOFF  /)  624-0566 

GO  TO  129  624-0567 

128  LSTAR  ■  V  /  AP(NSTAR)  624-0568 

IF  (LSTAR. GT.  10. 9143*( ALPHAN/RIL  )**2 )**. 4198)  GO  TO  130  624-0569 

WRITE  (3,105)  LSTAR  624-0570 

105  FORMAT  ( 1H0/46X.13HL*  WARNING  AT  F6.2/)  624-0571 

129  TAIL  *  -l.  624-0572 

130  WRITE  (3,104)  TIME,  P(L),  PIN),  P(N*l),  FG,  WIN),  624-0573 

1  (WEBP(I),  I  *  1 »N, I  PR ) ,  LOOPC,  P(NSTAR),  PTSTAR,  ISP,  SUMF,  SUMW,  624-0574 
1  DEFL (NSTAR)  624-0575 

104  FORMAT  ( 1H0F 5 .2 , 3F 8 . 1 , F 8 .0 ,  F8 . 3 ,9X  ,  10F7. 3, 1 6/2X , 3F8 . 1 , F8 .0 , 2F8 .3  )  624-0576 


(LIST.NE.O)  WRITE  (3,100)  LIST,  TIME,  FG,  ADEFL,  P(L) 

RETAIN  NOMINAL-CASE  PARAMETERS 
CALCULATE  SLOPE  MATRIX  TERMS 

(STEP.LT.O.)  GO  TO  90 
( TIME! .LT . TPART )  GO  TO  90 
=  IS  +  1 


624-0577 

624-0578 

624-0579 

624-0580 

624-0581 

624-0582 

624-0583 

624-0584 


A-  18 


o  o  o 
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IF  (STtP.GT.O.)  GO  TO  88 


624-0585 


COUNTP  =  COUNT  P  +  l. 

CALL  TAB  (COUNT P,  0..  T PART  ,  11) 
mi,  IS)  «  T  I  ME 
Hll.IS-Fl)  *  1.E6 
H  ( 2  ,  I  S  )  =■  PIC) 

H( 3, IS  )  «  FG 
H(4, IS  )  *  R ( N ) 

HI  5, IS  )  *  RINXRAV ) 

GO  TU  SO 

88  TPART  -  H I  I  ,  I  S  + 1  ) 

HI6.IS)  =  PIL)  /  HI2,  IS)  -  1. 

HI7, IS  )  *  FG  /  HI  3, IS)  -  1. 

H I  8 , I S  )  =  RIN)  /  H ( 4 ,  I S  )  -  1. 

HIS,  IS)  =  RINXRAV)  /  HI  5,  IS )  -  1. 

SO  IF  I T  IMEl.GT.TSTOP)  GO  TO  109 
IF  (PIL  )  .LT.PAMB+2. )  GO  TO  10S 
FGL  *  FG 
WNL  «  MIN) 

ML  =  M(N*1) 

VOLL  =  VOL 
PL  *  PIN+l) 

XIN)  *  XIN)  -  DELL 
IF  (MN.LE.O.)  GU  TO  91 

XDOT  =  I EXPMNT  *  (XTERN1  -  MN  *  I XTERH3  ♦  PLO  *  XDO )  «■  PDOT  ) 
1  ♦MOOT)  /  KO 

XOO  *=  XD  /  KD 
PLO  *  PIL) 

91  IF  (L.GT.l)  GO  TO  95 
XHEAO  =  XU) 

IF  I T IME.GT.TRISE )  OP  =  APT  /  APINSTAR) 

IF  ( T  I  ME.LE. TR1SE )  OP  *  IFIX-M.)  /  FIX 
DO  93  I  =L  , NS  TAR 

93  PII)  =  PI  I )  *  DP 

94  IF  (L.GT.l)  GO  TO  95 

CALL  TAB  (TIME,  0.,  OELT,  5) 

MF2  =  0. 

GO  TO  96 

95  L  *  L  -  1 

PIL)  =  PIL+l)  *  (FIX  ♦  1.)  /  FIX 
DELT  =*  I X (L  +  l )  -  X  ( L  )  )  /  VFZ 

96  LOOPC  «  -2 
FIX  *  FIX  ♦  1. 

APT  =  APINSTAR) 

TIME  *  TIME  ♦  DELT 
GO  TO  21 

OUTPUT  AND  ADVANCE  OF  PARAMETRIC  STUDIES 

109  TEMP  =  SUMF  /  FILL 

MR  I T  E  (3,108)  TEMP,  FILL 


624-0586 
624-0587 
624-0588 
624-0589 
624-0590 
624-0591 
624-0592 
624-0593 
624-0594 
624-0595 
624-0596 
624-0597 
624-0598 
624-0599 
624-0600 
624-0601 
624-0602 
624-0603 
624-0604 
624-0605 
624-0606 
624-0607 
624-0608 
624-0609 
624-0610 
624-0611 
624-0612 
624-0613 
624-0614 
624-0615 
624-0616 
624-0617 
624-0618 
624-0619 
624-0620 
624-0621 
624-0622 
624-0623 
624-0624 
624- 0625 
624-0626 
624-0627 
624-0628 
624-0629 
624-0630 
624-061 1 
624-0632 
624-0633 
624-0634 
624-0635 
624-0636 
624-0637 
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108  FORMAT  <12hO  STUP  ISP  «  F7.2»2HH,  BASEO  ON  LOADED  WEIGHT  OE  F7.1, 
1  AH  LbS  /  ) 

IE  ISTEP)  1,110,111 

110  WRITE  (3,1061  (H(l,l),  H( 2 ,  I  ) •  H(3,ll,  HI4,l»,  HIS, I),  1*1, IS) 

106  EOKHAT  ( 19H1NUM  INAL  PARAMETERS  /  IHQ4X ,  4HTIME  8X,  6hP  HEAD  6X, 
l  2HEG  10X,  2  EM  *  9X,  5HRIN4I  /  (5F12.3)) 

GO  TO  112 

111  Z(NPAR)  *  l <  NEAR )  /  1.1 

WRITE  13,107)  NPAR ,  (H(l, II,  H(6,i),  H( 7, I ) ,  H(8,I),  H(9,I), 

1  1*1, IS) 

107  FORMAT  I 39H1SENS I T I V  1  T  Y  DERIVATIVES  FOR  PARAMETER  13  /  1H04X, 

1  4hT  I  ME  8X,  6HP  HEAD  6X,  2HFG  10X,  2HR*  9X,  5HRIN4)  /  (5F12.3)) 

112  STEP  =  STEP  ♦  1. 

CALL  TAB  (STEP,  0.,  TEMP,  12) 

N  PAR  =  TEMP  l.E-6 
IF  (NPAR.EO.O)  GO  TO  1 
Z ( NPAR )  *  1.1  *  l ( NPAR  ) 

GO  TO  4 
END 


624-0638 

624-0639 

624-0640 

624-0641 

624-0642 

624-0643 

624-0644 

624-0645 

624-0646 

624-0647 

624-0648 

624-0649 

624-0650 

624-0651 

624-0652 

624-0653 

624-0654 

624-0655 

624-0656 
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C 

c 

c 


c 

c 


c 

c 

c 

c 


c 

c 

c 


SUBROUTINE  TAB  IN 

INPUT  LOGIC  EUR  VARIOUS  TABULAR  OATA 

U  I  MENS  I  ON  GEOMET!  3000)  ,  ML  1ST  (102),  LOGICI25) 

COMMON  /COM1/COM( 102  I  /C0M4/MA S TER <  3000)  /C0M5/SYSTEMI 125 ) 
EUUIVALENCE  (GEOMET, MASTER)  ,  (MUST, COM),  (  LOG  IC,  SYST  EM  ( 1 0 1  >  > , 
1  <J,ML1ST(101>),  (M, ML IST( 102 ) ) 


J  =  MAXO  <J,  l) 

2  ML  *  J 

CALL  OATA  IN  (0) 

I  Type  ■  LOGIC ( A ) 

IF  ( ITYPE.LT.O)  RETURN 
M  *  M  ♦  1 
MUST  (M)  *  ML 
IE  (M.LE.50)  GO  TO  1 
WRITE  (3,12) 

12  FORMAT  ( 1H010X, 35H  LIMIT  OF  50  TABLE  INPUTS  EXCEEDED  /) 
CALL  EXIT 
1  MX  *  LOGIC!  1) 

MUST  I MX  +  50 )  *  M 
NX  <  MAXO  (LOGIC ( 6) ,  II 
NY  *  MAXO  (LOGIC ( 5 )  ,  1) 

N  LOG  *  LOGIC! 3) 

MASTER(ML*3)  «  I  TYPE 
MASTER  (ML «-13l  *  NY 
MASTER! ML+2 )  *  N  LOG 

READ  -Y-  ARGUMENTS  (LIMIT  TO  19) 

MAKE  ALTERATIONS  FOR  NON-RECTANGUL AR  TABLES 

J  *  ML  ♦  1A 
INY  *  1 
JUMP  «  0 

if  ( i type )  a, a, 3 

3  INY  =  NY 
JUMP  *  1 

A  00  5  1*1,1  NY 

MASTER(J)  *  MAXO  (LOGICt 1+5) ,  1) 

5  J  1  J  ♦  1 

MAST ER ( ML  +  7 )  *  J 

00  7  1*1, NY 

GEOMET ( J )  *  SYSTEM!  1 ) 

IF  (N  LOG.GT.O)  GEOMET ( J  )  =  ALOG  ( GEOMET ( J ) ) 

7  J  *  J  ♦  1 
INOX  ■  J 

MASTER(ML+6)  *  INDX 
GEOMET ( ML+1 1 )  =  SYSTEM!  1 > 

REAO  -X-  ARGUMENTS  (LIMIT  TO  501 


624-0657 

62A-0658 

62A-0659 

62A-0660 

624-0661 

624-0662 

624-0663 

62A-066A 

624-0*65 

624-0666 

624-0667 

62A-0668 

624-0669 

62A-0670 

624-0671 

624-0672 

624-0673 

624-0674 

624-0675 

624-0676 

624-0677 

624-0678 

624-0679 

624-0680 

624-0681 

624-0682 

624-0683 

624-0684 

624-0685 

624-0686 

624-0687 

624-0688 

624-0689 

624-0690 

624-0691 

624-0692 

624-0693 

624-0694 

624-0695 

624-0696 

624-0697 

624-0698 

624-0699 

624-0700 

624-0701 

624-0702 

624-0703 

624-0704 

624-0705 

624-0706 

624-0707 

624-0708 

624-0709 
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Ji  ■  ML  ♦  14 

624-0710 

00  8  K-1.INY 

624-0711 

IF  (NY.GT.lt  CALL  DATA  IN 

il) 

624-0712 

NX  *  MASTERUl) 

624-0713 

00  6  1*1, NX 

624-0714 

IF  (N  L0G.GT.0I  SYSTEM!  1)  - 

ALOG 

•SYSTEM! 1 ) > 

624-0715 

GEOMET(J)  *  SYSTEM!  1  ) 

624-0716 

6 

J  ■  J  ♦  1 

624-0717 

8 

Jl  -  Jl  ♦  JUMP 

624-0718 

GEOMET!ML+10I  *  GEOMET ! INOX ) 

624-0719 

624-0720 

KEAO  TABLE  VALUES  -Z- 

624-0721 

624-0722 

1 NOZ  ■  J 

624-0723 

MASTERIML+B)  *  1N0Z 

624-0724 

Jl  *  ML  ♦  14 

624-0725 

00  10  K=1,NY 

624-0726 

CALL  DATA  IN  (1) 

624-0727 

NX  «  MASTER ( J 1 ) 

624-0728 

00  9  1*1, NX 

624-0729 

IF  (N  LOG.GT.O)  SYSTEM!  I)  - 

ALOG 

(SYSTEM! It) 

624-0730 

GEOMETIJI  =  SYSTEM! It 

624-0731 

9 

J  *  J  ♦  1 

624-0732 

10 

Jl  *  Jl  ♦  JUMP 

624-0733 

GEOMET ( ML*12 )  *  GEOMET! INOZt 

624-0734 

GO  TO  2 

624-0735 

ENO 

624-0736 

o  O  <~l  <1  O  O 
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C 

C 

C 

C 

c 

c 

c 


c 


c 

c 

c 


c 


SUBROUTINE  TAB  1X1,  Y1 ,  Zl,  MX) 


624-0737 


624-0738 

TABLE  LOOK-UP  ROUTINE  FOR  LINEAR  INTERPOLATION  OF  624-0739 

ONE-OI MENS  I ONAL  Z»F(X),  OR  Z-FIY),  TWO  DIMENSIONAL  Z*F,X,Y),  624-0740 
ANO  TWO-DIMENSIONAL  NON-RECTANGUL AR  TABLES  WITH  COLUMNS  OF  624-0741 
VARYING  NUMBERS  OF  ROWS  624-0742 

PERFORM  LOG-LOG  INTERPOLATION  IF  CODED  (MASTER! 3)  OR  -(ML*2)>  624-0743 

624-0744 


DIMENSION  GEOMET ( 3000 ) t  ML  1ST ( 102  )  ,  A(2),  JUNM32) 

COMMON  /COMl/COMI  102)  /C0M4/MASTER  (  3000)  /C0M5/SYSTEM(  125 ) 
EQUIVALENCE  (GEOMET, MASTER) ,  (MLIST.COM),  ( JUNK, SYSTEM! 90) ) 
EQUIVALENCE  ( I  TYPE, JUNKI  4) ) ,  ( JXl , JUNK ( 5 > >  ,  I JV1, JUNK ( 6 ) ) , 

1  (  I NOX  * JUNK (7)1,  ( INDY,JUNK(8) ),  ( INDZ, JUNK(9) ) , 

2  (XS,JUNK(11 )),  (  YS , JUNK (121),  ( ZS , JUNK (13)), 

3  (  NY , JUNK (14)1,  (NX1 , JUNK ( 15 ) ) 

M  ■  ML  1ST  (MX*50> 

IF  (Ml  9,7,10 

CODE  TO  ALLOW  ONLY  ONE  ERROR  MESSAGE 


624-0745 

624-0746 

624-0747 

624-0748 

624-0749 

624-0750 

624-0751 

624-0752 

624-0753 

624-0754 

624-0755 

624-0756 


624-0757 


7 

WRITE  (3,8)  MX,  M 

624-0758 

8 

FORMAT  I40H0  REFERENCE  TO  UNCORRELATED  TABLE  NUMBER  214 

/  ) 

624-0759 

MLIST  ( MX+50 )  ■  -l 

624-0760 

9 

21  -  0. 

624-0761 

RETURN 

624-0762 

624-0763 

10 

X  *  XI 

624-0764 

Y  *  Y1 

624-0765 

ML  *  MLIST(M) 

624-0766 

IF  (MASTERIML+2 ) . EQ.O )  GO  TO  25 

624-0767 

X  *  ALOG  (X) 

624-0708 

Y  «  ALOG  (Y) 

624—0769 

624-0770 

MOVE  INDICES  FROM  PREVIOUS  ENTRY  TO  WORKING  STORAGE 

-JUNK- 

624-0771 

624-0772 

25 

LOW  *  0 

624-0773 

J  «  ML 

624-0774 

DELY  ■  0. 

624-0775 

DO  27  1*1,30 

624-0776 

JUNKI I)  *  MASTER (J ) 

624-0777 

27 

J  •  J  ♦  1 

624-0778 

NX  *  NX1 

624-0779 

624-0780 

DETERMINE  COLUMN  AND  Y  INTERPOLATION  RATIO  -DEL  Y- 

624-0781 

SHORT  CIRCUIT  IF  INPUT  ARGUMENTS  ARE  UNCHANGED  FROM 

PRIOR 

SET624-0782 

624-0783 

IF  (ABSIXS  -  X)  ♦  ABSIYS  -  Yl)  67,67,31 

624-0784 

31 

DYl  *  Y  -  GEOMET (INOY) 

624-0785 

IF  I 0Y1 )  38,50,33 

624-0786 

32 

OYl  *  DY2 

624-0787 

33 

IF  (JY1  ♦  1  -  NY)  34,50,50 

624-0788 

34 

0Y2  *  Y  -  GEOMET! INDY+1 ) 

624-0789 

A-23 


sO 


oooo  o  r'  r~.  o  o 
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IF  ( DY2  I  40,35,35 

SEQUENCE  FORWARD  ONE  COLUMN 

35  INDY  «  INDY  ♦  1 
IF  (ITYPE)  37,37,36 

36  NX  *  JUNK ( JYI +  1 5  I 
INUX  =  INOX  ♦  NX 

37  INOZ  «  I NDZ  ♦  NX 
JYI  =  JYI  +  I 
IF  (LON)  32,32,50 

SEQUENCE  BACKWARDS  ONE  COLUMN 

38  IF  (JYI)  50,50,39 

39  INDY  ■=  INDY  -  1 
JYI  =  JYI  -  1 
IF  ( ITYPE)  42,42,41 

41  NX  =  J  UNK ( J  Y 1  + 1 5  I 
INUX  =  I NDX  -  NX 

42  INDZ  =  INDZ  -  NX 
DY2  =  DYl 

DYl  =  Y  -  GEUMET ( INDY) 

IF  (DYl)  38,50,40 

40  DELY  *  DYl  /  (DYl  -  DY2I 

CONTINUE  TO  FIND  X  INTERPOLATION  RATIO  AND  Z  VALUE  (ONE  OR 
TWO  VALUES  AS  REQUIRED) 

50  DELX  >=  0. 

C  COUNT  BACK  TO  THE  END  OF  THE  COLUMN  ( NON-RECTANGULAR » 

IF  (ITYPE)  51,51,48 

48  NX  =  JUNK  (  J  Y  1  1  5 ) 

I  =  JXl  ♦  1  -  NX 
IF  (I)  51,51,49 

49  JXl  =  JXl  -  I 

I  NDX  a  I  NDX  -  I 
INOZ  *  INDZ  -  I 
C 

51  DX1  ■■  X  -  GEOMET  (  INOX ) 

IF  (DX1)  55,58,52 

52  IF  (JXl  ♦  1  -  NX)  53,58,58 

53  0X2  =  X  -  GEOMET ( INDX  +  1  ) 

IF  ( 0X2 )  57,54,54 

C  COUNT  FORWARD  ONE  ROW 

54  JXl  =  JXl  ♦  1 
INUX  »  I NDX  ♦  l 
INOZ  «  INDZ  ♦  1 
0X1  =  0X2 

GO  TO  52 

55  IF  (JXl)  58,56,56 

C  COUNT  BACKWARDS  ONE  ROW 

56  JXl  =  JXl  -  1 
INOX  =  I  NDX  -  1 
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IN02  -  INOZ  -  1 
0X2  -  0X1 

DX1  ■  X  -  GEOMET ( INDX ) 

IF  I  0X1  )  $5,58,57 

57  OELX  *  OXt  /  (0X1  -  0X2) 

A  (LOW* l )  GEOHET  (  INOZ  )  ♦  OELX  *  I GEOMET  ( 1 NDZ+ 1 )  -  GEOMET  ( INDZ  ) ) 
GO  TO  60 

58  A  I L  OW+1 )  *  GEOMET! INOZ) 

RESTORE  STORED  INDICES  AND  SET-UP  SECOND  X  INTERPOLATION 
SECONO  X  VALUE  AND  A  Y  INTERPOLATION  IS  NEEDED  ONLY  IF 
-OELY-  HAS  BEEN  SET  POSITIVE 

60  IF  (LOW)  61,61,65 

61  LOW  *  1 

J  *  ML  «■  4 
DO  62  1*5,9 

MASTER  I J )  *  JUNK ( 1 ) 

62  J  *  J  ♦  l 

IF  IDELY)  63,63,35 

63  GEOMET IML+12)  *  A(l) 

GO  TO  66 

COMPLETE  TWO-DIMENSIONAL  INTERPOLATION 

SAVE  VALUES  X,  Y,  AND  Z  FOR  SHORT  CIRCUIT  ON  RE-ENTRY 

65  GEOME  T ( ML  +  12 )  *  A(l)  ♦  OELY  *  <A(2)  -  All)) 

66  GEOMET (ML+10)  *  X 
GEOMET (ML  +  l 1 )  *  Y 

67  ANS  *  GEOMETIML+12 ) 

Z1  *  ANS 

IF  (MASTERIML+2).  NE.O)  Z1  =  EXP  (ANS) 

RETURN 

END 
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SUBROUTINE  OATA  IN  (Nl 

REAO  RANDOM-ENTRY  DATA  TO  GSTA8LISH  PROBLEM  OEF  iNITION 

STACK  EACH  GROUP  WITH  INTEGERS  FIRST,  THEN  FLOATING  OATA 
CONTROL  KEY  -N-  CALLS  FOR  OATA  ONLY  IF  (♦),  BOTH  IF  10),  OR 
INTEGERS  ONLY  IF  (-1 

READ  UP  TO  6  PAIRS  (LOCATION  12,  OATA  (14,6X1  OR  FIO)  PER  CARO 

ANY  FIELO  MAY  BE  BLANK  -  NUMBER  FIRST  LOCATION  IN  LISTS 

(->  LOCATION  OR  BLANK  CARD  TERMINATES  SEQUENCE 

OIMENS ION  LIT),  K(6),  E(6),  L0GICC25),  COMENT (7) 

COMMON  /COM 5/ SYSTEM! 125  » 

EQUIVALENCE  ( LOGIC ,SYST£M(101 )) 

IF  (N)  1,1,11 
00  2  1*1,25 

LOGIC(I)  *  0 

REAO  (1,3)  (L( I ),K(  I ), 1*1,4),  L(5),  (COMENTf I ) , 1-1 ,71 
FORMAT  (4(12,14,  6X),  12,  7A4) 

WRITE  (3,231  (COMENT(I), 1-1,7) 

FORMAT  ( /6H  LOGIC  10X,  7A4 I 

J  *  0 

K1  *  1 

00  9  1*1,4 

IF  (  L(  I ) )  5,6,5 

J  -  L(  II 

IF  (J)  10,8,8 

IF  (Jl  7,9,7 

J  *  J  ♦  1 

IF  (K(II)  8,9,8 

LOGIC! J)  *  Kill 

K1  *  I 

CONTINUE 

WRITE  (3,4)  (L(I>,K(II,I~1,K1),  L(K1*1) 

FORMAT  (6X,6( 18,15, 7X)I 
IF  (Jl  30,30,29 
IF  ( L( 5 ) )  30,22,22 
IF  (Nl  21,11,11 

00  12  1*1,50 

SYSTEM! I)  *  0. 

READ  (1,13)  (L( t), E(  11,1*1,61,  L ( 7) 

FORMAT  (6( 12, FIO. 5)  ,  121 
J  *  0 
K1  *  1 

00  19  1-1,6 

IF  ( L( 1 1  I  15,16,15 

J  *  L(  1 1 

IF  (Jl  20,18,18 
IF  (Jl  17,19,17 
J  *  J  ♦  1 

IF  ( E( 1 1 1  18,19,18 
SYSTEM! J I  *  E (  I  ) 

K1  *  I 
CONTINUE 

WRITE  (3,141  (L(II,E(II,I*1,K1I,  L(K1>1) 

FORMAT  (6H  OATA  6(I8,F12.6),  141 
IF  (Jl  21,21,24 
IF  ( L(7) )  21,32,32 

RETURN 

END 
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